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1,  Introduction 


The  androgen  receptor  (AR)  plays  an  essential  role  in  prostate  cancer,  and  suppression  of  its 
signaling  with  androgen  deprivation  therapy  (ADT)  has  been  the  mainstay  of  treatment  for 
metastatic  hormone  sensitive  prostate  cancer  for  more  than  70  years.  We  now  appreciate  that  the 
AR  remains  a  critical  factor  activated  by  residual  androgens  in  castration  resistant  prostate  cancer 
(CRPC)  patients.  The  residual  prostatic  androgen  in  part  or  largely  derived  from  uptake  of 
circulating  adrenal  androgens.  DHEA  (dehydroepiandrosterone)  and  its  sulfated  form,  DHEAS 
are  secreted  in  large  amounts  by  the  adrenal  cortex  and  are  precursors  for  the  production  of  T 
(testosterone)  andDHT  (dihydrotestosterone),  the  most  potent  androgens,  in  the  tumor 
microenvironment.  Genetic  variations  at  two  SLCO  (Solute  Carrier  Organ  Anion  Transporter 
Family)  genes  (3  SNPs  for  SLC02B1  and  2  SNPs  for  SLC01B3)  may  alter  androgen  uptake 
affecting  prostate  cancer  outcomes.  SLC02B1/SLC01B3  exhibit  higher  expression  levels  in 
CRPC  tumors  than  in  normal  or  early  stage  prostate  cancers,  suggesting  these  genes  are 
important  in  maintaining  the  CRPC  phenotype.  Statins  and  estrogens  have  been  shown  to  enter 
cells  using  the  SLC02B1/SLC01B3  transporters.  Statin  usage  appears  to  have  a  significant 
impact  on  the  progression  to  advanced  prostate  cancer.  Estrogens  have  an  antitumor  effect  in 
CRPC.  We  hypothesized  that  statins  and  estrogens  might  potentially  block  the  uptake  of 
androgen  precursors  from  the  tumor  microenvironment  via  interactions  with 
SLC02B1/SLC01B3.  Therefore,  elucidating  the  role  and  inhibiting  the  transport  activity  of 
SLC02B1/SLC01B3  in  CRPC  may  ultimately  lead  to  improved  outcomes  of  prostate  cancer. 

We  intend  to  demonstrate  that  statins/estrogens  compete  for  DHEAS/T  uptake  through  the 
SLC02B1/SLC01B3  transporters  of  in  prostate  cancer.  Combining  ADT  and  statins  or  inhibitors 
of  SLC02B1/SLC01B3  may  significantly  enhance  the  efficacy  of  ADT  and  the  overall  survival. 
Our  study  will  reveal  mechanisms  mediating  the  development  of  CRPC;  enhance  our 
understanding  of  molecular  mechanisms  underlying  the  impact  of  statins  or  estrogens  on  prostate 
cancer;  and  suggest  new  therapeutic  targets  potentially  improving  the  efficacy  of  ADT. 

2,  Keywords 

Prostate  cancer,  androgen  deprivation  therapy  (ADT),  time  to  progression  (TTP),  overall  survival 
(OS,)  SLC02B1,  SNP,  DHEAS,  androgen,  statins. 

3,  Accomplishments 

►  What  were  the  major  goals  of  the  project? 

Aim  1:  We  will  determine  the  impact  of  statins  and  estrogens  on  the  uptake  of  DHEAS  orT 
through  SLC02B1  and  SLC01B3  and  their  effects  on  prostate  cancer  cell  lines. 

Aim  2:  We  will  determine  the  impact  of  expression  levels  of  SLC02B1  and  SLC01B3  on 
prostate  tumor  progression  to  CRPC  in  the  presence  or  absence  of  statins  or  estrogens  in  animal 
models. 
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Aim  3:  We  will  determine  the  association  between  prostate  cancer  survival  and  the  use  of  statins 
and  tumor  expression  levels  of  SLC02B1  and  SLC01B3  in  men  with  different  genotypes  of 
SLC02B1  and  SLC01B3  in  the  Health  Professionals  Follow-up  Study  and  ADT  prostate  cancer 
cohorts. 

Aim  4:  We  will  identify  and  optimize  small  molecule  inhibitors  of  SLC02B1  and  SLC01B3  and 
test  these  inhibitors  in  cell  based  and  animal  models  of  prostate  cancer. 

►  What  was  accomplished  under  these  goals? 

A.  During  the  past  few  years,  we  studied  statin  use  at  the  time  of  initiation  of  androgen 
deprivation  therapy  (ADT)  and  time  to  progression  (TTP)  in  patients  with  hormone-sensitive 
prostate  cancer  (Mentioned  in  Aim  1)  and  examined  the  association  of  SLC02B1  genotypes  with 
TTP  and  overall  survival  (OS)  in  patients  receiving  ADT  for  prostate  cancer  (Mentioned  in  Aim 
3). 

Since  statins  utilize  SLC02B1  to  enter  cells,  we  hypothesized  that  they  may  compete  with 
DHEAS  uptake  and  that  statin  use  might  prolong  TTP  while  on  ADT  (Aim  1).  To  determine  if 
statins  interfere  with  DHEAS  uptake,  we  performed  in  vitro  studies  using  prostate  cancer  cell 
lines.  Next,  we  queried  our  institutional  clinical  database  to  assess  for  an  association  between 
statin  use  and  TTP  on  ADT. 

Using  our  IRB  approved  institutional  clinical  database  (Prostate  CRIS),  we  identified  1265 
patients  with  hormone-sensitive  PC  who  had  been  treated  with  ADT  (with  or  without  an 
antiandrogen)  from  January  1996  until  November  2013.  Patients  were  excluded  if  they  had 
insufficient  follow-up  data  on  PSA  after  ADT  administration  (n=131)  or  if  statin  use  status  was 
unknown  (N=208),  which  left  926  patients  for  this  analysis.  Clinicodemo graphic  data  was 
captured  from  the  CRIS  database.  The  electronic  medical  record  was  reviewed  for  dates  of 
initiation  and  progression  on  ADT  and  whether  patients  were  on  statins  at  the  time  of  ADT 
initiation.  Progression  was  defined  as  a  minimum  of  3  PSA  rises.  Date  of  progression  was 
defined  as  date  of  first  PSA  rise  (nadir  +  >0.02  ng/ml)  or  radiologic  progression.  Patient  and 
disease  characteristics  were  summarized  as  frequencies  or  the  median  and  range  of  values. 
Characteristics  were  compared  between  statin  users  and  non-users  using  chi-square  and 
Wilcoxon  rank  sum  tests.  The  primary  outcome  variable  was  TTP  on  ADT,  defined  as  the 
duration  of  time  from  ADT  initiation  to  the  date  of  disease  progression  or  censored  at  the  date  of 
last  follow-up  visit  in  patients  who  were  progression-free.  The  association  between  statin  use  and 
TTP  on  ADT  was  estimated  from  multivariable  Cox  regression  to  estimate  hazard  ratios  (HR) 
and  95%  confidence  intervals  (Cl),  adjusting  for  pre-defined  prognostic  factors:  biopsy  Gleason 
score,  primary  therapy  type,  use  of  prior  ADT  in  conjunction  with  local  therapy,  metastatic  status, 
and  PSA  at  ADT  initiation. 

Patient  and  disease  characteristics  at  diagnosis  and  at  ADT  initiation  are  detailed  in  Table  1  of 
Appendices  Paper  1  (page  41).  Statin  users  tended  to  have  a  lower  median  PSA  both  at 
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diagnosis  (9.1  vs.  11.8)  and  at  ADT  initiation  (10.3  vs.  12.5).  Duration  from  diagnosis  to  ADT 
initiation  was  longer  in  statin  patients  (3.9  vs.  2.3  years).  Users  were  more  likely  to  have  lower 
stage  disease  (56%  vs.  44%  T1  disease)  and  less  likely  to  have  de  novo  metastases  (11%  vs.  18%) 
or  nodal  involvement  (5%  vs.  10%)  at  diagnosis.  Statin  users  were  more  likely  to  have 
undergone  local  therapy  or  to  have  received  ADT  as  part  of  local  therapy  (33%  vs.  26%). 

Patients  on  statins  were  less  likely  to  have  metastases  at  ADT  initiation  (53%  versus  63%)  (p 
values  <0.05). 

At  the  time  of  data  capture,  70%  (n=644) 
of  patients  had  progressed  on  ADT  by 
PSA.  Median  follow-up  was  5.8  years 
(range:  0.1-15.9).  Median  TTP  on  ADT 
for  all  patients  irrespective  of  statin  use 
was  20.3  months  (95%  Cl:  17.5,  23.6). 

Statin  users  at  ADT  initiation  had  a 
significantly  longer  median  TTP  on 
ADT  (27.5  vs.  17.4  months,  p=0.0005, 

Fig  1A).  Since  some  but  not  all  of  the 
baseline  characteristics  favored  more 
indolent  disease  in  the  statin  users,  we 
adjusted  for  these  potential  imbalances. 

Importantly,  the  association  remained 
statistically  significant  after  adjusting 
for  pre-defined  prognostic  clinical 
factors  including  biopsy  Gleason  score, 
type  of  primary  therapy,  use  of  prior 
ADT  in  conjunction  with  localized  therapy,  metastatic  status  and  PSA  at  initiation  of  ADT 
(adjusted  HR=0.83,  95%  Cl:  0.69,0.99,  Table  2  of  Appendices  Paper  1  (page  44)).  When 
stratified  by  year  of  ADT  initiation  using  5-year  increments,  the  association  between  statin  use 
and  TTP  on  ADT  remained  significant  in  the  multivariable  model  (adjusted  HR=0.83,  95%  Cl: 
0.69,1.00).  Moreover,  the  association  between  statin  use  and  TTP  was  observed  regardless  of 
whether  patients  had  radiographic  evidence  of  metastatic  disease  compared  to  biochemical 
relapse  only  at  ADT  initiation  (HR=0.79,  95%  CI:0.58,  1.07  for  M0;  HR=0.84,  95% 

0:0.67,1.06  for  Ml,  p  for  interaction=0.72).  (Fig.  IB,  1C) 

In  our  ADT  cohort  of  926  patients,  283  (31%)  were  taking  a  statin  at  ADT  initiation.  After  a 
median  follow-up  of  5.8  years,  644  patients  (70%)  had  progressed  on  ADT.  Median  TTP  on 
ADT  was  20.3  months  (95%  Cl:  18,24).  Men  on  statins  had  a  longer  median  TTP  on  ADT 
compared  to  non-users  (27.5  vs.  17.4  months,  p=0.0005).  The  association  remained  statistically 
significant  after  adjusting  for  pre-defined  prognostic  factors  [adjusted  HR=0.83  (p=0.039)].  The 
positive  statin  effect  was  observed  for  both  patients  with  and  without  metastases.  Therefore, 
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Figure  1:  Kaplan  Meier  analysis  of  TTP  on  ADT  according 
to  statin  use  in  all  patients.  Kaplan  Meier  plot  of  TTP  on 
ADT  according  to  statin  use  (A)  In  all  patients  (B)  In  patients 
with  no  visible  radiographic  metastasis  but  biochemical  (PSA) 
failure  (MO)  and  (C)  in  patients  with  metastasis  (Ml)  at  ADT 
initiation. 
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statin  use  at  the  time  of  ADT  initiation  was  associated  with  a  significantly  longer  TTP  on  ADT 
even  after  adjusting  for  known 
prognostic  factors. 

In  vitro,  we  demonstrated  that  statins 
block  DHEAS  uptake  by  competitively 
binding  to  SLC02B1.  We  examined  the 
effect  of  four  different  statins 
(atorvastatin,  fluvastatin,  pravastatin  and 
simvastatin)  on  DHEAS  uptake  in  the 
androgen  dependent  LNCaP  and  the 
partially  androgen  dependent  22RV 1  PC 
cell  lines  (Fig.  2A).  DHEAS  uptake  in 
PC  cell  lines  was  concentration  and  time 
dependent  (Fig.  2B,  2C).  When 
incubated  with  DHEAS  at  a 
physiological  concentration  (2.5  pM)  for 
60  min,  the  22RV 1  cell  line,  which  has  a 
relatively  high  level  of  SLC02B 1 
expression,  displayed  the  most  active 
DHEAS  uptake  of  more  than  300 
pmol/mg  compared  to  ~60  pmol/mg 
protein  for  LNCaP  (Fig.  2B).  100  pM 
atorvastatin  significantly  decreased 
DHEAS  influx  by  -50%  in  both  cell  lines 
when  cells  were  incubated  with  2.5  pM 
DHEAS.  Among  the  four  statins  we 
studied,  pravastatin  had  the  most 
significant  inhibitory  effect  on  DHEAS 
uptake  in  both  cell  lines,  while  a  more 
prominent  effect  of  simvastatin  was 
shown  for  22RV1  cells  than  that  of  LNCaP  cells  (Fig.  2B).  However,  10  or  even  100  pM 
atorvastatin  or  simvastatin  was  insufficient  to  inhibit  DHEAS  uptake  in  LNCaP,  which  has  a 
relatively  low  level  of  SLC02B1  expression  when  the  concentration  of  100  pM  DHEAS  was 
used.  These  results  suggest,  not  surprisingly,  that  different  statins  compete  with  DHEAS  for  the 
same  transporter  but  with  varying  efficiency  and  that  this  effect  is  cell  line  dependent. 

To  determine  the  dependence  of  DHEAS  uptake  on  SLC02B1,  we  constructed  inducible 
SLC02B 1  deficient  stable  22RV 1  cell  lines  using  the  lentiviral  derived  tetracycline  inducible 
shRNA  knock  down  system.  We  were  unable  to  establish  an  inducible  SLC02B 1  deficient  stable 
cell  line  in  LNCaP,  thus,  we  used  a  transient-inducible  knock  down  of  SLC02B 1  LNCaP  cells  in 
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Figure  2:  Statins  inhibit  DHEAS  uptake  by  SLC02B1 
in  PC  cells.  A.  Relative  mRNA  levels  in  PC  cells  before 
and  after  SLC02B1  is  knocked  down.  B.  DHEAS  uptake 
in  PC  cells  with  2.5  pM  DHEAS  and  different 
concentrations  of  statins  when  incubated  for  60  minutes. 
Statistical  analysis  was  performed  by  comparing  each 
condition  with  the  DHEAS  2.5  pM  and  no  statin  state 
(NS)  except  when  indicated.  C.  Depicts  DHEAS  uptake 
in  PC  cells  before  (scrambled  shRNA  (scr))  and  after 
(shRNA  2B1)  SLC02B1  is  knocked  down  when 
incubated  with  2.5  pM  DHEAS  and  100  pM  atorvastatin 
(ATO)  for  10  and  60  minutes.  Statistical  analysis  was 
performed  by  comparing  each  condition  with  scrambled 
shRNA  (scr)  after  10  min  with  DHEAS  except  when 
indicated. 
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this  study  (shRNA-SLC02Bl).  After  successfully  knocking  down  SLC02B1  (Fig.  2A),  DHEAS 
uptake  was  substantially  decreased  to  -50%  in  22RV1  and  -70%  in  LNCaP  of  that  observed  in 
control  cells  (cells  transfected  with  scrambled  shRNA)  (Fig.  2C).  These  results  indicate  that 
SLC02B1  plays  an  essential  role  in  DHEAS  import  into  PC  cells  (Fig.  2C).  More  importantly, 
knocking  down  SLC02B1  abolished  atorvastatin’s  inhibition  of  DHEAS  uptake;  further 
suggesting  that  atorvastatin  competes  with  DHEAS  for  binding  to  their  transporter,  SLC02B 1 . 

To  further  support  our  hypothesis,  we  wanted  to  demonstrate  that  inhibition  of  cell  growth  was 
SLC02B 1  and  DHEAS  dependent.  Thus,  we  examined  the  impact  of  atorvastatin,  the  most 
commonly  used  statin  clinically,  and  DHEAS  on  tumor  proliferation  before  and  after  SLC02B 1 
was  knocked  down  using  the  WST-1 
assay. 

2.5  pM  (or  a  higher  concentration) 
atorvastatin  is  known  to  inhibit  LNCaP 
cell  proliferation  and  induces 
autophagy.  Since  atorvastatin 
concentrations  in  patient  serum  range 
from  5-270  nM,  we  chose  200  nM  of 
atorvastatin  for  our  cell  proliferation 
assay.  80  nM  DHEAS  significantly 
increased  cell  proliferation  in  LNCaP 
and  22RV 1  lines  that  were  maintained 
in  androgen-depleted  medium.  Cell 
numbers  nearly  doubled  by  day  6  for 
LNCaP  (Fig.  3A).  At  day  6,  DHEAS 
induced  a  -  6-fold  increase  in  LNCaP 
cell  number,  compared  to  an  only  -3 
fold  increase  in  the  absence  of  DHEAS 
(p=  0.0003).  However,  treatment  with 
atorvastatin  (200  nM)  inhibited  this 
DHEAS  induced  cell  proliferation. 

Consistent  with  this  finding,  knocking  down  SLC02B1  abolished  DHEAS  induced  cell 
proliferation  in  LNCaP  and  22RV1  (Fig.  3B).  Furthermore,  treatment  with  atorvastatin  did  not 
significantly  inhibit  the  growth  of  SLC02B1  knocked-down  cells.  Taken  together,  our  data 
demonstrate  that  atorvastatin  can  efficiently  block  SLC02B  l-mediated  DHEAS  uptake  and 
DHEAS  induced  cell  growth  in  androgen  dependent  PC  cell  lines. 

Our  in  vitro  findings  that  statins  competitively  reduce  DHEAS  uptake  and  thus,  effectively 
decrease  the  available  intratumoral  androgen  pool,  affords  a  plausible  mechanism  to  support  the 
clinical  observation  of  prolonged  TTP  in  statin  users.  This  study  was  published  in  JAMA 
oncology  (Please  see  the  Appendices  Paper  1,  page  19). 


Figure  3.  Atorvastatin  (ATO)  decreases  DHEAS  induced 
PC  cell  proliferation  through  SLC02B1.  Atorvastatin 
(ATO)  decreases  DHEAS  induced  PC  cell  proliferation  by 
competing  for  influx  through  SLC02B1.  A.  WST-1  assay  of 
cell  proliferation  in  PC  cell  lines.  B.  WST-1  assay  of  cell 
proliferation  in  PC  cell  lines  after  SLC02Blis  knocked  down. 
For  proliferation  assays,  cells  were  maintained  in  androgen 
depleted  medium  with  doxycycline  for  2  days  followed  by  the 
addition  of  80  nM  DHEAS  and  200  nM  ATO  to  the  culture 
medium  for  the  indicated  time.  Relative  cell  numbers  were 
calculated  as  percentages  of  the  cell  numbers  at  day  0  (100%). 
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In  addition,  we  sought  to  validate  the  association  of  three 
previously  demonstrated  germline  variants  in  SLC02B1 
with  TTP  on  ADT  and  to  evaluate  if  these  genetic 
variants  were  associated  with  OS  for  prostate  cancer 
(partial  of  Aim  3).  The  three  single  nucleotide 
polymorphisms  (SNPs)  were  genotyped  in  an 
independent  validation  cohort  of  616  patients  with  PC 
treated  with  ADT  at  our  institution  from  1996  to  2013. 

Multivariable  Cox  regression  adjusting  for  known 
prognostic  factors  estimated  the  association  of  genetic 
variants  with  TTP  on  ADT  and  OS.  The  expression  of 
SLC02B 1  was  examined  in  prostatectomy  samples.  The 
impact  of  SLC02B 1  expression  level  on  uptake  of 
DHEAS  was  evaluated  in  cell  lines. 

In  this  study,  the  clinicodemographic,  disease,  and  past 
treatment  characteristics  of  the  initial  (N=478)  and 
validation  ADT  cohort  (n=616)  are  presented  in 
Supplementary  Table  1  of  Appendices  Paper  2 
(page76).  The  two  cohorts  had  similar  disease 
characteristics  at  diagnosis  and  at  ADT  initiation.  In  both 
cohorts,  approximately  70%  of  patients  had  received  a 
local  therapy  (radical  prostatectomy  or  radiation  therapy) 
and  nearly  60%  had  metastases  at  the  time  of  ADT 
initiation.  The  validation  cohort  had  a  lower  PSA  at 
diagnosis  (median  9.9  vs.  14  ng/ml),  probably  reflecting 
the  more  frequent  use  of  PSA  screening  leading  to  an 
earlier  diagnosis  of  PC  in  those  patients.  The  validation  cohort  more  frequently  received 
intermittent  ADT  (35.9%  vs  18.8%). 

We  first  validated  the  Association  of  SLC02B 1  SNPs  with  TTP  during  ADT  in  the  validation 
cohort  (N=616).  At  the  time  of  data  retrieval,  66%  (n=408)  of  patients  in  the  validation  cohort 
had  progressed  on  ADT.  The  median  TTP  on  ADT  was  20.9  months  (95%  Cl:  18.0,  24.0)  and 
the  median  follow-up  time  was  4.2  years  (range:  0.1, 16.3). 

The  exonic  SNP,  rs  12422 149,  in  SLC02B1  was  associated  with  TTP  on  ADT  in  univariable 
analysis  (median  TTP:  27.2  (95%  Cl:  18.9,48.9)  months  for  AA/AG  and  20.0  (16.5,23.0)  months 
for  GG,  P  =  .019)  (Fig  4  and  Table  1  of  Appendices  Paper  2  (page  68)).  In  multivariable 
analyses  adjusting  for  clinical  factors,  the  association  between  the  SNP  and  TTP  remained 
significant  (HR=1.31,  95%  Cl:  1.00,  1.72,  P=.049).  Thus,  the  association  of  rsl2422149  with 
TTP  on  ADT  observed  in  the  initial  cohort  was  validated  such  that  PC  patients  carrying  the 
major  GG  genotype  exhibited  a  shorter  TTP  on  ADT.  While  there  was  a  trend  toward  an 
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rs12422149 
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-  GG 


association  between 
rs  1077858  with  TTP  on 
ADT,  the  two  intronic  SNPs 
were  not  validated  (Fig  4B, 

4C). 

We  then  explored  the 
correlation  between  the 
SLC02B1  SNPs  and  TTP  in 
the  combined  initial  and 
validation  patients  and 
stratified  by  metastatic 
disease  status  at  time  of  ADT  initiation,  both  of  which  had 
not  been  previously  evaluated.  In  the  combined  cohort 
(n=1094),  74%  of  patients  (n=81 1)  had  progressed  on 
ADT.  Median  TTP  on  ADT  was  18.9  months  (95%  Cl: 
16.5,  21.1).  For  the  exonic  rsl2422149,  we  again  found  a 
significant  association  with  TTP  on  ADT  in  the  combined 
cohort  (adjusted  HR=1.33,  95%CI:  1.10-1.60,  P=  .003). 
When  stratified  by  metastatic  disease  status,  the 
association  remained  in  the  Ml  population  (adjusted 
HR=1.65,  95%:  1.28-2.12)  but  was  not  seen  in  the  M0 
population  (adjusted  HR=0.96,  95%  Cl:  0.72-1.27)  ( P 
interaction=  .006,  Fig  5  and  Table  2  of  Appendices  Paper 
2  (page  69)).  The  association  between  the  intronic  SNP 
rs  1077858  and  TTP  was  of  borderline  significance 
( P=  .075;  P=  .062  if  AA  and  AG  were  combined).  No 
association  was  found  for  rs  1789693  and  TTP  on  ADT  in 
the  combined  cohort.  Results  were  similarly  negative  in  the 
M0  and  Ml  populations  {P  interaction  >0.5).  The 
associations  with  TTP  on  ADT  were  similar  among 
patients  with  and  without  prior  hormone  treatments  for  the 
exonic  SNP  rs  12422 149.  However,  for  the  intronic  SNP 
rsl077858,  a  significant  association  was  observed  only  in 
patients  without  prior  hormone  treatment  (adjusted 
HR=1.32  (1.06,  1.65)),  but  not  in  patients  with  prior 
hormone  treatment  (adjusted  HR=0.87  (0.55,  1.37)). 

Correlations  with  OS  had  not  been  explored  in  our  prior 
work  and  thus  we  evaluated  this  endpoint  in  the  combined 
cohort.  Nearly  half  (49%,  n=537/1094)  of  the  patients  had 
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Figure  5.  Kaplan-Meier  curves  of  time  to  progression  during  ADT  in  all 
cohorts  (initial  plus  validation  cohort)  according  to  rsl2422149  genotypes 
as  stratified  by  metastatic  status  (A)  M0  versus  (B)  Ml  at  ADT  initiation. 
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Fig  6.  Kaplan-Meier  curves  of  overall 
survival  from  ADT  initiation  in  all 
patients  (initial  plus  validation 
cohort),  according  to  SLC02B1 
genotypes  (A)  rsl2422149,  (B) 

rsl789693,  (C)  rsl077858. 
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died  at  the  time  of  data  collection.  The  median  OS  from  ADT  initiation  was  6.5  years  (95%  Cl: 
6.0,  7.0)  and  the  median  follow-up  time  was  6.5  years  (range  0.1,  16.3  years). 

There  was  no  statistically  significant  association  with  OS  for  the  exonic  SNP  rs  12422 149  in 
either  univariable  or  multivariable  analyses  (Fig  6A  and  Table  3  of  Appendices  Paper  2  (page 
70)).  In  univariable  analysis,  there  was  no  association  between  rsl789693  and  OS  from  ADT 
initiation  ( P=  .184,  Bonferroni 
adjusted  P=0.552),  but  patients 
carrying  the  minor  allele  (AT  or 
TT)  had  longer  OS  in 
multivariable  analysis  (HR=0.81 
and  0.78,  respectively,  P=  .044) 

(Fig  6B  and  Table  3 
Appendices  Paper  2  (page  70)). 

Importantly,  we  found  that 
patients  having  the  minor 
genotype  GG  in  the  intronic  SNP 
rs  1077858  had  a  shorter  OS  from 
ADT  initiation  in  both 
univariable  ( P=  .009,  Bonferroni 
adjusted  P= 0.027)  and 
multivariable  analysis  (adjusted 
HR=1.35(1.07,  1.71),  P=  .012) 

(Fig  6C  and  Table  3 
Appendices  Paper  2  (page  70)). 

The  median  OS  was  decreased 
from  6.7  (95%  Cl:  6.2,  7.2)  to 
5.2  (95%  Cl:  4.3,  6.8)  years  with 
the  AA/AG  vs.  GG  rs  1077858 
genotypes  respectively. 

Given  the  fact  that  the  intronic 
SNP  rs  1077858  was  associated 
with  OS,  we  studied  the  possible 
mechanism  by  which  the  SNP 
might  affect  OS  from  ADT 
initiation.  Using  80  normal 
prostate  tissue  samples  available 
in  our  tissue  repository  of 
patients  known  to  have  PC,  we 
analyzed  SLC02B 1  expression 


Figure  7.  Impact  of  the  intronic  SNP  rsl077858  on  the  expression 
of  SLC02B1  which  affects  the  sulfated  dehydroepiandrosterone 
(DHEAS)  uptake.  (A)  Relative  SLC02B1  mRNA  levels  compared  to 
GAPDH  in  normal  prostate  tissue  in  primary  prostate  cancer  patients. 
Lines  in  the  figure  represent  the  median  with  interquartile  range.  (B) 
Relative  SLC02B1  mRNA  levels  compared  to  GAPDH  in  22RV1 
cells  with  different  alleles  after  CRISPR  reconstitution.  Lines  in  the 
figure  represent  the  median  with  interquartile  range.  (C)  SLC02B1 
mRNA  levels  in  mock  transfected  LNCaP  cells  and  LNCaP  cells 
transfected  with  different  SLC02B1  variants.  Values  represent  the 
fold  differences  relative  to  those  in  mock  transfected  cells,  which  were 
set  as  1.0. WT,  LNCaP  mock  transfected  with  either  siRNA  control  or 
pCMV6-XL4  vector;  siSLC02Bl,  LNCaP  transiently  transfected  with 
siRNA  targeting  SLC02B1  (-p  represents  the  SMARTpool  siRNA 
from  GE  Healthcare;  -1,  -2  and  -3  represent  three  unique  27mer 
siRNAs,  respectively.);  oeSLC02Bl,  LNCaP  transiently  transfected 
with  pCMV6-XL4-SLC02Bl.  Statistical  analysis  (unpaired  t-test) 
was  performed  by  comparing  each  condition  with  the  WT  condition. 
(D)  DHEAS  uptake  in  LNCaP  cells.  Cells  were  incubated  with  2.5 
pM  DHEAS  for  10  minutes  after  cells  have  been  transfected  with 
siRNA  or  pCMV6-XL4-SLC02Bl  for  3  days.  All  experiments  were 
repeated  in  triplex.  Statistical  analysis  (unpaired  t-test)  was  performed 
by  comparing  each  condition  with  the  WT  condition. 
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(Fig  7A).  Results  demonstrated  that  the  different  rsl077858  SNPs  were  associated  with 
significantly  different  SLC02B1  expression  levels  (Ptrend=  .019)  with  the  highest  SLC02B1 
expression  in  patients  carrying  the  minor  GG  genotype.  To  prove  that  different  rsl077858  SNPs 
impact  the  expression  of  SLC02B1,  we  reconstituted  the  rsl077858  SNP  variants  using  CRISPR 
in  22RV 1  cells  to  convert  the  original  A  allele  to  the  G  allele  (Supplementary  data).  The 
association  between  different  rs  1077858  SNPs  and  SLC02B1  expression  levels  was  significant 
(P trend=  .0073)  and,  once  again,  the  GG  genotype  in  22RV1  cells  corresponded  to  the  highest 
SLC02B1  expression  level  (Fig  7B).  These  results  suggest  that  the  association  of  rs  1077858  on 
OS  was  possibly  through  changes  in  SLC02B1  expression  levels  in  these  patients.  As  a  result, 
alleles  corresponding  to  higher  SLC02B1  expression  may  allow  greater  DHEAS  uptake  and 
thereby  induce  shorter  OS  in  patients.  To  evaluate  this  hypothesis,  we  examined  DHEAS  uptake 
in  a  LNCaP  PC  cell  line  by  altering  the  expression  level  of  SLC02B1  to  mimic  different 
rsl077858  allele  conditions.  After  successfully  knocking  down  or  overexpressing  SLC02B1  (Fig 
7C),  we  found  that  DHEAS  uptake  activity  was  dependent  on  the  expression  level  of  SLC02B1 
and  that  higher  expression  of  SLC02B1  transported  more  DHEAS  into  the  cells  (Fig 
7D). Therefore,  the  association  of  the  SNP  rsl077858  with  OS  may  be  in  part  due  to  differential 
SLC02B1  expression  and  consequent  higher  uptake  of  DHEAS  and  subsequent  resistance  to 
ADT,  which  in  turn  might  contribute  to  shortened  survival.  This  study  was  published  in  JCO 
(Please  see  the  Appendices  Paper  2  (page  48)). 

B.  Our  study  of  the  function  of  SLCO  transporters  suggested  that  SLC02B1  may  potentially 
\import  abiraterone  acetate  (AA).  We  evaluated  the  effects  of  AA  and  statins  alone  and  in 
combination  in  PC  cell  lines  and  studied  the  impact  of  statin  use  on  AA  treatment  duration  in 
patients  with  CRPC. 

The  key  finding  was  that  AA’s  effect  on  PC  cell  proliferation  is  SLC02B1  dependent.  When 
SLC02B1  is  knocked  down  (shSLC02Bl)  in  both  22RVland  LNCAP  cells,  AA’s  inhibition  of  cell 
growth  is  attenuated  (Figure  8).  These  findings  suggest  that  AA  efficacy  depends  on  SLC02B 1 
transport  into  cells.  To  assess  for  competition  between  AA  and  statins  for  uptake,  they  evaluated 
the  effect  of  AA  on  22RV1  cell  line  proliferation  and  assessed  whether  it  was  dependent  on 
SLC02B1.  Atorvastatin  attenuated  AA’s  inhibition  of  the  22RV1  cell  proliferation  and  was 
SLC02B1  dependent.  Consistent  with  this  observation,  data  from  the  Mostaghel  Lab 
demonstrated  decreased  ABI  levels  in  OATP2B1 -expressing  LNCaP  cells  in  the  presence  of 
fluvastatin  (Figure  8C).  Across  the  different  statins,  fluvastatin  and  simvastatin  seemed  to 
attenuate  AA’s  inhibition  of  cancer  growth  the  most  with  pravastatin  being  less  inhibitory 
(Figure  8  E,F). 
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Figure  8.  (A)  AA  efficacy  is  dose  and  SLC02B 1  dependent.  When  SLC02B1  is  knocked  down 
(sh)  in  22RV1  cells,  AA’s  inhibition  of  cell  growth  is  attenuated.  (B)  The  greatest  degree  of  AA 
inhibition  occurs  in  the  absence  of  statin  but  in  the  presence  of  intact  SLC02B1.  In  the  presence  of 
statin,  inhibition  of  cell  proliferation  by  AA  (indicated  by  ABI)  is  decreased  and  is  concentration 
dependent.  When  SLC02B1  is  knocked  down,  AA’s  inhibitory  effects  on  cell  growth  are 
attenuated.  (C)  Impact  of  statin  on  AA  levels  in  VC  (vector  control)  vs  SLC02B1 -expressing  cells. 
Cells  treated  with  or  without  fluvastatin  (FLU)  lOOuM  for  10  minutes  prior  to  ABI  0.5uM  for  4 
hours.  (D)  ABI  (AA)  levels  in  LNCaP  cells  overexpressing  the  indicated  SLCO  transporters  vs. 
vector  control  (VC).  Cells  were  treated  with  ABI  0.5uM  for  4  hours.  (E,  F)  LNCaP  cells  were 
treated  with  different  concentration  of  ABI  (no  abi  treatment,  2.5  pM,  5.0  pM)  and  with  or  without 
various  statins  (ATO,  atorvastatin;  FLU;  PRA,  pravastatin;  SIM,  simvastatin)  at  10  nM  and  200 
nM.  Relative  cell  numbers  were  examined  by  WST-1  assay  at  day  6.  The  value  for  each  column 
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represents  the  percentage  of  cell  numbers  comparing  to  the  no  ABI  treatment  condition  within  each 
statin  treatment  group.  Fluvastatin  induced  the  greatest  attenuation  of  ABI  effects  at  the  higher 
concentration  (200  nM)  whereas  pravastatin  showed  the  least  inhibition  of  ABI  effects,  especially 
when  ABI  was  used  near  the  physiological  concentration  (2.5  pM). 

Thus,  we  hypothesized  that  statins  may  compete  with  AA  for  influx  by  SLC02B 1,  which  could 
negatively  impact  drug  efficacy.  We  queried  our  IRB  approved  institutional  clinical  database  for 
patients  who  had  received  AA  for  mCRPC.  Data  on  statin  use,  including  on-therapy  dates,  were 
retrospectively  collected.  Duration  on  AA  was  used  as  a  surrogate  for  TTP. 

Of  the  224  patients  treated  with  AA  between  2008  and  2015,(GU,  ASCO,  2016)  nearly  half  were 
statin  users  (41%).  There  was  a  surprising  trend  to  longer  AA  duration  in  statin  users:  14.2  vs.  9.2 
months  (HR  0.79,  95%  Cl:  0.57-1.09,  p=0.14,  Figure  9).  The  difference  was  even  more  pronounced 
in  patients  who  were  docetaxel  or  enzalutamide  naive  (21.3  vs.  14.8  mo,  HR  0.75,  Cl:  0.49,  1.15). 
PSA  at  AA  start  was  similar  between  the  non-users  and  users  (20.3  vs.  16.9  ng/dL),  and  differences 
were  not  thought  to  be  due  to  baseline  imbalances  in  the  patient  populations. 


Without  statin  use 

With  statin  use  -  -  -  - 


Figure  9.  Trend  to  longer 
duration  of  AA  use  in  statin  users 
(top)  compared  to  non-users 
(bottom)):  14.2  vs.  9.2  months 
(HR  0.79) 


In  summary,  contrary  to  our  initial  hypothesis,  there  was  a  trend  to  longer  AA  duration  in 
statin  users.  Given  our  prior  work  demonstrating  a  median  10  month  prolonged  TTP  on  ADT  in 
men  concurrently  taking  statins  (similar  adjusted  HR  0.83)  and  in  vitro  findings  of  inhibition  of 
DHEAS  uptake  by  statins,  we  postulate  that  there  could  be  an  additive  effect  with  statins  blocking 
DHEAS  uptake  and  AA  inhibiting  androgen  synthesis. 


►  What  opportunities  for  training  and  professional  development  has  the  project  provided? 


The  project  has  provided  opportunities  for  post-doctoral  training  and  for  the  collaboration 
between  physician  and  scientists.  The  success  of  this  project  has  helped  the  participated  post- 


11 


doctoral  fellow  moving  onto  his  next  carrier  level  and  the  participated  physician  developing  new 
considerations  in  prostate  cancer  therapies. 

►  How  were  the  results  disseminated  to  communities  of  interest? 

We  found  that  statin  use  at  the  time  of  ADT  initiation  was  associated  with  a  significantly  longer 
TTP  on  ADT,  most  likely  due  to  its  competition  on  the  DHEAS  uptake.  The  beneficial  effect  of 
statin  usage  has  broadly  drawn  attentions  in  prostate  cancer  field. 

►  What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

•  Determine  the  impact  of  estrogens  on  the  uptake  of  DHEAS  or  T  through  SLC02B 1  and 
SLC01B3  and  their  effects  on  prostate  cancer  cell  lines. 

•  Determine  the  association  between  PCa  survival  and  the  use  of  statins  and  tumor 
expression  levels  of  SLC02B1  and  SLC01B3  in  men  with  different  genotypes  of  SLC02B1 
and  SLC01B3  in  the  Health  Professionals  Follow-up  Study  and  ADT  PCa  cohorts. 

•  Determine  the  impact  of  expression  levels  of  SLC02B1  and  SLC01B3  on  prostate  tumor 
progression  to  CRPC  in  the  presence  or  absence  of  statins  or  estrogens  in  animal  models. 

•  Demonstrate  the  impact  of  statins  on  SLCO  transporters-mediated  transport  of  DHEAS  and 
AA  in  prostate  tumors  and  the  extent  to  which  these  effects  influence  prostate  tumor  growth  in 
vivo  and  in  vitro  and  evaluate  the  addition  of  pravastatin  at  the  time  of  PSA  progression  on 
AA/pred  using  PSA  response  as  a  biologic  signal  of  antitumor  efficacy  in  patients  with  mCRPC 
in  a  prospective,  randomized  clinical  trial. 


4,  Impact 

►  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Almost  all  prostate  cancers  (PCa)  depend  on  androgens  for  growth  in  its  initial  stages.  We  can 
suppress  testicular  androgen  production  with  androgen  deprivation  therapy  (ADT),  which  is  the 
most  effective  and  widely  used  therapy  for  PCa  patients.  However,  residual  androgen  inpatients 
still  can  active  the  androgen  receptor  (AR)  and  lead  to  castration  resistant  PCa  (CRPC),  which  is 
usually  lethal.  DHEA  (dehydroepiandrosterone)  and  its  sulfated  form,  DHEAS  are  secreted  in 
large  amounts  by  the  adrenal  cortex  and  are  precursors  for  the  production  of  T  (testosterone)  and 
DHT  (dihydrotestosterone),  the  most  potent  androgens.  Therefore,  blocking  the  uptake  of  adrenal 
androgens  in  the  prostate  could  enhance  the  efficacy  of  ADT.  SLC02B1  and  SLC01B3  are  two 
transporter  genes  involved  in  DHEAS  and  T  uptake  in  the  prostate.  Evidence  exists  that 
SLC02B1  and  SLC01B3  are  expressed  at  higher  levels  in  advanced  PCa  and  the  expression 
levels  and  genotypes  of  these  two  genes  are  correlated  with  the  uptake  activity  of  DHEAS  and  T 
and  outcomes  of  PCa.  Since  statins  and  estrogens  have  been  shown  to  use  the  SLC02B1  and 
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SLC01B3  transporters  and  have  an  effect  on  the  advanced  prostate  cancer,  we  hypothesize  that 
statins  and  estrogens  might  work  by  blocking  the  uptake  of  androgens  via  interactions  with 
SLC02B1  and  SLC01B3.  Therefore,  illustrating  the  role  and  inhibiting  the  transport  activity  of 
SLC02B1  and  SLC01B3  in  CRPC  will  ultimately  lead  to  improve  outcomes  of  PCa.  In  this 
study,  we  will  determine  the  mechanism  if  statins/estrogens  compete  for  the  DHEAS/T  uptake 
through  the  SLC02B1  and  SLC01B3  transporters.  We  will  also  design  specific  inhibitors  of 
SLC02B1  and  SLC01B3.  A  combination  treatment  with  ADT  and  statins  or  inhibitors  of 
SLC02B1  and  SLC01B3  may  significantly  enhance  the  efficacy  of  ADT  and  outcomes  in  PCa 
patients. 

New  FDA  approved  agents  such  as  abiraterone  and  enzalutamide  successfully  target  the  androgen 
pathway.  But  resistance  in  patients  ultimately  develops  and  adaptive  mechanisms  permitting  the 
continued  stimulation  of  AR  by  persistent  ligand  occurs.  Novel  therapeutic  targets 
potentially  improving  the  efficacy  of  ADT  by  further  inhibition  of  residual  androgens  are  needed. 
Based  on  our  proposed  study,  we  plan  to  develop  SLC02B1  and  SLC01B3  genotypes  as 
biomarkers,  but  also  improve  the  outcomes  of  PCa  through  inhibition  of  SLC02B1  and 
SLC01B3.  Herein  we  will  validate  the  basis  of  the  putative  benefit  seen  with  statins  and 
develop  more  potent  inhibitors  of  this  important  adaptive  pathway.  If  experimental  results  in  our 
study  develop  as  anticipated,  we  estimate  that  within  3-5  years,  a  clinically  meaningful  outcome 
will  be  achieved. 

This  study  will  reveal  the  mechanisms  of  CRPC  development  and  enhance  our  understanding  of 
molecular  mechanisms  underlying  the  impact  of  statins  or  estrogens  on  PCa  progression.  It  will 
pave  a  new  road  to  improve  the  efficacy  of  ADT  by  focusing  on  new  therapeutic  targets  in  PCa 
patients. 

►  What  was  the  impact  on  other  disciplines? 

Nothing  to  report. 
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►  What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

►  What  was  the  impact  on  society  beyond  science  and  technology? 

It  draws  attentions  on  the  potential  effect  of  statin  use  on  tumor  progression. 

5,  Changes/Problems 

►  Changes  in  approach  and  reasons  for  change 

Nothing  to  report. 

►  Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  report. 

►  Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  report. 

►  Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents 

Nothing  to  report. 

►  Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  report. 

►  Significant  changes  in  use  or  care  of  vertebrate  animals. 

Nothing  to  report. 

►  Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  report. 

6,  Products 

►  Publications,  conference  papers,  and  presentations 

Harshman,  L.  C.,  Wang,  X.,  Nakabayashi,  M.,  Xie,  W.,  Valenca,  L.,  Werner,  L.,  Yu,  Y.,  Kantoff, 
A.  M.,  Sweeney,  C.  J.,  Mucci,  L.  A.,  Pomerantz,  M.,  Lee,  G.  S.,  and  Kantoff,  P.  W.  (2015)  Statin 
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Use  at  the  Time  of  Initiation  of  Androgen  Deprivation  Therapy  and  Time  to  Progression  in 
Patients  With  Hormone-Sensitive  Prostate  Cancer,  JAMA  oncology  1,  495-504. 

Wang,  X.,  Harshman,  L.C.,  Xie,  W.,  Nakabayashi,  M.,  Qu,  F.,  Pomerantz,  M.,  Lee,  G.S.,  and 
Kantoff,  P.W.  (2015)  Association  of  SLC02B1  Genotypes  with  Time  to  Progression  and  Overall 
Survival  in  Patients  Receiving  Androgen  Deprivation  Therapy  for  Prostate  Cancer,  Journal  of 
Clinical  Oncology.  Accepted. 

Philip  W.  Kantoff,  Xiaodong  Wang,  Wanling  Xie,  Mari  Nakabayashi,  Mark  Pomerantz,  Gwo- 
Shu  Mary  Lee.  Genetic  variants  of  the  organic  anion  transporter  SLC02B 1  as  pharmacogenomic 
determinants  of  response  to  androgen  deprivation  therapy  (ADT)  for  prostate  cancer. 
Presentation,  2015  Genitourinary  Cancers  Symposium. 

Lauren  Christine  Harshman,  Xiaodong  Wang,  Mari  Nakabayashi,  Wanling  Xie,  Loana  Bueno 
Valenca,  Lillian  Werner,  Yongjiang  Yu,  Aaron  Kantoff,  Christopher  Sweeney,  Lorelei  A.  Mucci, 
Mark  Pomerantz,  Gwo-Shu  Mary  Lee,  Philip  W.  Kantoff.  Statin  use  at  the  time  of  initiation  of 
androgen  deprivation  therapy  and  time  to  progression  in  patients  with  hormone-sensitive  prostate 
cancer.  Poster,  2015  ASCO  Annual  Meeting. 

Lauren  C.  Harshman,  Lillian  Werner,  Mari  Nakabayash,  Abhishek  Triphath,  Rana  McKay, 
Xiaodong  Wang,  Mark  Pomerantz,  Lorelei  A.  Mucci,  Mary-Ellen  Taplin,  Christopher  J. 
Sweeney,  Gwo-Shu  Mary  Lee,  Philip  W.  Kantoff.  The  impact  of  statin  use  on  abiraterone  acetate 
(AA)  treatment  duration  in  patients  with  castration-resistant  prostate  cancer  (CRPC).  2016  GU 
ASCO 

►  Website(s)  or  other  Internet  site(s) 

Statin  drugs  can  delay  prostate  cancer  progression  in  patients  receiving  androgen  deprivation 
therapy:  http://www.dana-farber.org/Newsroom/News-Releases/Statin-drugs-can-delay-prostate- 

cancer-progression-in-patients-receiving-androgen-deprivation-therapv-studv-shows.aspx 

Statins  Associated  with  Longer  Prostate  Cancer  Time  to  Progression  During  Androgen 
Deprivation  Therapy:  http://media.iamanetwork.com/news-itern/statins-associated-with-longer- 
prostate-cancer-time-to-progression-during-androgen-deprivation-therapv/ 

Statin  users  have  longer  prostate  cancer  time  to  progression: 
http://www.healthvlivingmagazine.us/Articles/9973/ 

Statin  drugs  can  delay  prostate  cancer  progression  in  patients  receiving  androgen  deprivation 
therapy,  study  shows:  http://www.sciencedailv.com/releases/2015/05/1505071140Q4.htm 

Concurrent  Statin  Use  in  ADT  Patients  Reduced  All  Cause,  Prostate  Cancer-Specific  Deaths: 
http://www.onclive.com/publications/urologists-in-cancer-care/2015/August-2Q15/concurrent- 
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statin-use-in-adt-patients-reduced-all-cause-prostate-cancer-specific- 

deaths#sthash.6BQd7wAk.dpuf 


►  Technologies  or  techniques 

Nothing  to  report. 


►  Inventions,  patent  applications,  and/or  licenses 

Nothing  to  report. 

►  Other  Products 

Nothing  to  report. 


7,  Participants  &  Other  Collaborating  Organizations 
►  What  individuals  have  worked  on  the  project? 


Name: 

Philip  Kantojf 

Project  Role: 

PI 

Researcher  Identifier: 

Nearest  person  month: 

12 

Contribution  to  Project: 

study  design,  preclinical  studies,  data  analysis,  data  intepretation 

Funding  Support: 

Dana-Farber  Prostate  Cancer  SPORE  P50CA090381 

Name: 

Xiaodong  Wang 

Project  Role: 

Key  personal 

Researcher  Identifier: 

Nearest  person  month: 

12 

Contribution  to  Project: 

preclinical  studies,  data  collection,  data  analysis,  data  intepretation 

Funding  Support: 

Name: 

Wanting  Xie 

Project  Role: 

Statistician 

Researcher  Identifier: 

Nearest  person  month: 

4 

Contribution  to  Project: 

statistical  analysis 

Funding  Support: 
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Name: 

Gwo-Shu  Mary  Lee 

Project  Role: 

Assistant  to  PI 

Researcher  Identifier: 

Nearest  person  month: 

12 

Contribution  to  Project: 

studies,  data  collection,  data  analysis,  data  intepretation 

Funding  Support: 

Name: 

Lauren  Harshman 

Project  Role: 

Researcher 

Researcher  Identifier: 

Nearest  person  month: 

6 

Contribution  to  Project: 

data  analysis,  data  interpretation 

Funding  Support: 
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►  Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key 
personnel  since  the  last  reporting  period? 

Yes. 

►  What  other  organizations  were  involved  as  partners? 

Nothing  to  report. 

8,  Special  Reporting  Requirements 

►  COLLABORATIVE  AWARDS 

Nothing  to  report. 

►  QUAD  CHARTS 
Nothing  to  report. 
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9,  Appendices 
Paper  1: 

Statin  use  at  the  time  of  initiation  of  androgen  deprivation  therapy  and  time  to  progression 
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Abstract: 


Importance:  Statin  use  has  been  associated  with  improved  prostate  cancer  outcomes.  DHEAS 
is  a  precursor  of  testosterone  and  a  substrate  for  SLC02B1,  an  organic  anionic  transporter.  We 
previously  demonstrated  that  genetic  variants  of  SLC02B1  correlated  with  time  to  progression 
(TTP)  on  androgen  deprivation  therapy  (ADT). 

Objective:  Since  statins  utilize  SLC02B1  to  enter  cells,  we  hypothesized  that  they  may  compete 
with  DHEAS  uptake  and  that  statin  use  might  prolong  TTP  while  on  ADT. 

Design:  In  vitro  studies  using  prostate  cancer  cell  lines.  A  retrospective  analysis  of  prospectively 
collected  data  on  926  patients,  who  had  received  ADT  from  January  1996  until  November  2013. 
Setting:  Academic,  comprehensive  cancer  center. 

Participants:  Patients  who  had  received  ADT  for  biochemical  or  metastatic  recurrence  or  de 
novo  metastatic  prostate  cancer. 

Main  Outcomes  and  Measures:  To  determine  if  statins  interfere  with  DHEAS  uptake,  we 
performed  in  vitro  studies  using  prostate  cancer  cell  lines.  Next,  we  queried  our  institutional 
clinical  database  to  assess  for  an  association  between  statin  use  and  TTP  on  ADT.  TTP  was 
estimated  using  multivariable  Cox  regression  and  adjusted  for  known  prognostic  factors. 

Results:  In  vitro,  we  demonstrated  that  statins  block  DHEAS  uptake  by  competitively  binding  to 
SLC02B1.  In  our  ADT  cohort  of  926  patients,  283  (31%)  were  taking  a  statin  at  ADT  initiation. 
After  a  median  follow-up  of  5.8  years,  644  patients  (70%)  had  progressed  on  ADT.  Median  TTP 
on  ADT  was  20.3  months  (95%  Cl:  18,24).  Men  on  statins  had  a  longer  median  TTP  on  ADT 
compared  to  non-users  (27.5  vs.  17.4  months,  p=0.0005).  The  association  remained  statistically 
significant  after  adjusting  for  pre-defined  prognostic  factors  [adjusted  HR=0.83  (p=0.039)].  The 
positive  statin  effect  was  observed  for  both  patients  with  and  without  metastases. 
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Conclusions  and  Relevance:  Statin  use  at  the  time  of  ADT  initiation  was  associated  with  a 


significantly  longer  TTP  on  ADT  even  after  adjusting  for  known  prognostic  factors.  Our  in  vitro 
findings  that  statins  competitively  reduce  DHEAS  uptake  and  thus,  effectively  decrease  the 
available  intratumoral  androgen  pool,  affords  a  plausible  mechanism  to  support  the  clinical 
observation  of  prolonged  TTP  in  statin  users. 
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Background 


The  organic  anionic  transporter,  SLC02B1  enables  a  variety  of  anticancer  compounds  and 
hormones  to  enter  cells.1  Amongst  its  substrates  is  the  abundant  adrenal  androgen 
dehydroepiandrosterone  sulfate  (DHEAS),  which  is  a  precursor  to  more  potent  androgens,  such 
as  dihydroxytestosterone  (DHT),  which  binds  to  the  androgen  receptor  (AR)  in  normal  and 
cancer  cells.  In  prostate  cancer  (PC),  expression  of  SLC02B1  increases  upon  progression  from 
hormone  sensitive  to  castration-resistant  disease.'  Our  group  and  others  have  previously 
demonstrated  that  genetic  variants  in  SLC02B1  are  associated  with  the  durability  of  response  to 
androgen  deprivation  therapy  (ADT),  due  to  varied  efficiency  of  androgen  influx  into  cells.3’4 

Interestingly,  statins  are  also  substrates  of  SLC02B1.  Past  work  has  generally  shown  an  inverse 

c  o 

association  between  statin  use  and  incidence  of  PC  as  well  as  improved  clinical  outcomes.  ’ 

Little  is  known  about  the  impact  of  statin  use  and  the  durability  of  response  to  ADT,  which  is  the 
cornerstone  of  treatment  for  metastatic  hormone  sensitive  PC.9  Given  the  fact  that  both  DHEAS 
and  statins  are  substrates  for  SLC02B1,  we  first  sought  to  determine  if  there  was  any  interaction 
between  statins  and  DHEAS  influx  by  SLC02B1  in  PC  cell  lines.  Then,  using  our  prospectively 
collected  institutional  database,  we  evaluated  the  association  between  statin  use  and  time  to 
progression  (TTP)  among  PC  patients  on  ADT. 

Methods 

Cell  lines  and  reagents 

The  hormone  sensitive  PC  cell  lines,  LNCaP  and  22RV1  were  used  in  the  study.  LNCaP  and 
22RV1  cells  were  maintained  in  RPMI  1640  and  supplemented  with  10%  FBS  and  antibiotics. 
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For  the  cell  proliferation  studies,  all  PC  cells  were  cultured  in  Phenol-Red  free  RPMI 1640/10% 
charcoal- stripped  FBS.  293T  cells  were  obtained  from  the  American  Type  Culture  Collection 
(ATCC)  and  maintained  in  Dulbecco's  Modified  Eagle's  Medium  (DMEM)  with  10%  FBS  and 
antibiotics.  All  cell  lines  were  regularly  screened  for  mycoplasma  (Sigma  Venor  GeM 
Mycoplasma  Detection  Kit).  DHEAS  was  obtained  from  BioVender  R&D  Products.  Atorvastatin 
was  obtained  from  Santa  Cruz  Biotechnology,  and  fluvastatin,  pravastatin  and  simvastatin  were 
purchased  from  Selleckchem. 

shRNA 

We  constructed  an  inducible  shSLC02Bl  RNA  expressing  plasmid  using  the  “all-in-one” 
pLKO-Tet-On  lentiviral  vector  (see  supplementary  eMethods).  Lentiviruses  were  packaged  using 
293T  cells.  The  stable  inducible  shSLC02Bl -expressing  cell  lines  were  established  by  selection 
in  puromycin.  The  efficiency  of  knocking  down  SLC02B1  expression  was  assayed  after 
induction  with  1  ng/ml  of  doxycycline  for  48  hours.  The  LNCaP  and  22RV1  cell  lines  were 
transfected  with  scrambled  shRNA,  and  these  were  used  as  a  negative  controls. 

Quantitative  RT-PCR 

100  ngs  total  RNA  was  extracted  from  each  PC  cell  line  and  analyzed  by  RT-PCR  (see 
supplementary  eMethods).  All  RT-PCR  experiments  were  performed  in  triplicate. 

DHEAS  uptake  assay 

Cells  were  harvested  in  the  PBS  buffer  and  triple  washed  with  incubation  buffer  (140  mM  NaCl, 
5  mM  KC1,  ImM  KH2PO4,  1.2  mM  MgSCE,  1.5  mM  CaCh,  5  mM  D-glucose,  and  12.5  mM 
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HEPES,  pH  7.7).  Aliquoted  cells  were  incubated  with  each  statin  or  DMSO  (control)  in 
incubation  buffer  at  37°C  for  10  minutes,  and  then  treated  with  different  DHEAS  concentrations. 
At  defined  time  points,  the  cells  were  triple  washed  with  cold  incubation  buffer  to  stop  DHEAS 
uptake.  The  cells  were  then  lysed  using  1%  Triton  X-100  solution  in  1  x  PBS  on  ice  for  30 
minutes.  The  total  protein  concentration  was  measured  in  the  cleared  cell  lysate  supernatant  by 
protein  BCA  (bicinchoninic  acid)  assay  (Thermo  Scientific).  The  quantity  of  intracellular 
DHEAS  was  determined  using  the  DHEAS  ELISA  kit  (BioVendor,  Candler,  NC)  and  adjusted  to 
the  protein  concentration  of  the  cell  lysate. 

Cell  proliferation  assay 

Cell  proliferation  was  determined  using  the  WST-1  assay  (Roche,  Indianapolis,  IN).  Briefly, 
control  and  shSLC02Bl-lentivirus  infected  cells  were  cultured  in  96-well  plates  in  the  presence 
of  doxycycline  at  a  confluence  of  -10%  in  androgen-depleted  medium  for  2  days  followed  by 
treatment  with  atorvastatin  and/or  DHEAS.  Cell  proliferation  assays  were  carried  out  on 
different  days  after  treatment.  Each  experiment  was  performed  in  triplicate. 

Clinical  cohort  study 

Using  our  IRB  approved  institutional  clinical  database  (Prostate  CRIS),10  we  identified  1265 
patients  with  hormone-sensitive  PC  who  had  been  treated  with  ADT  (with  or  without  an 
antiandrogen)  from  January  1996  until  November  2013.  Patients  were  excluded  if  they  had 
insufficient  follow-up  data  on  PSA  after  ADT  administration  (n=131)  or  if  statin  use  status  was 
unknown  (N=208),  which  left  926  patients  for  this  analysis. 
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Clinicodemographic  data  was  captured  from  the  CRIS  database.  The  electronic  medical  record 
was  reviewed  for  dates  of  initiation  and  progression  on  ADT  and  whether  patients  were  on 
statins  at  the  time  of  ADT  initiation.  Progression  was  defined  as  a  minimum  of  3  PSA  rises.  Date 
of  progression  was  defined  as  date  of  first  PSA  rise  (nadir  +  >0.02  ng/ml)  or  radiologic 
progression. 

Statistical  analysis 

Patient  and  disease  characteristics  were  summarized  as  frequencies  or  the  median  and  range  of 
values.  Characteristics  were  compared  between  statin  users  and  non-users  using  chi-square  and 
Wilcoxon  rank  sum  tests.  The  primary  outcome  variable  was  TTP  on  ADT,  defined  as  the 
duration  of  time  from  ADT  initiation  to  the  date  of  disease  progression  or  censored  at  the  date  of 
last  follow-up  visit  in  patients  who  were  progression-free.  The  association  between  statin  use  and 
TTP  on  ADT  was  estimated  from  multivariable  Cox  regression  to  estimate  hazard  ratios  (HR) 
and  95%  confidence  intervals  (Cl),  adjusting  for  pre-defined  prognostic  factors:  biopsy  Gleason 
score,  primary  therapy  type,  use  of  prior  ADT  in  conjunction  with  local  therapy,  metastatic  status, 
and  PSA  at  ADT  initiation.11'12 

For  in  vitro  studies,  data  was  represented  as  means  ±  S.D.  of  at  least  3  biological  repeats. 
Comparison  between  two  independent  groups  (or  cell  lines)  was  performed  by  a  2-tailed  t  test.  P 
<  0.05  was  considered  statistically  significant  for  all  analyses. 
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Results 


Inhibition  of  DHEAS  uptake  by  atorvastatin  is  SLC02B1  dependent 

We  examined  the  effect  of  four  different  statins  (atorvastatin,  fluvastatin,  pravastatin  and 
simvastatin)  on  DHEAS  uptake  in  the  androgen  dependent  LNCaP  and  the  partially  androgen 
dependent  22RV1  PC  cell  lines  (Fig.  1A).  DHEAS  uptake  in  PC  cell  lines  was  concentration  and 
time  dependent  (Fig.  IB,  1C,  and  eFig.  1).  When  incubated  with  DHEAS  at  a  physiological 
concentration  (2.5  pM)  for  60  min,  the  22RV1  cell  line,  which  has  a  relatively  high  level  of 
SFC02B1  expression,  displayed  the  most  active  DHEAS  uptake  of  more  than  300  pmol/mg 
compared  to  ~60  pmol/mg  protein  for  LNCaP  (Fig.  IB).  100  pM  atorvastatin  significantly 
decreased  DHEAS  influx  by  -50%  in  both  cell  lines  when  cells  were  incubated  with  2.5  pM 
DHEAS.  Among  the  four  statins  we  studied,  pravastatin  had  the  most  significant  inhibitory 
effect  on  DHEAS  uptake  in  both  cell  lines,  while  a  more  prominent  effect  of  simvastatin  was 
shown  for  22RV1  cells  than  that  of  LNCaP  cells  (Fig.  IB).  However,  10  or  even  100  pM 
atorvastatin  or  simvastatin  was  insufficient  to  inhibit  DHEAS  uptake  in  LNCaP,  which  has  a 
relatively  low  level  of  SLC02B1  expression  when  the  concentration  of  100  pM  DHEAS  was 
used  (eFig.l).  These  results  suggest,  not  surprisingly,  that  different  statins  compete  with 
DHEAS  for  the  same  transporter  but  with  varying  efficiency  and  that  this  effect  is  cell  line 
dependent. 

To  determine  the  dependence  of  DHEAS  uptake  on  SLC02B1,  we  constructed  inducible 
SLC02B1  deficient  stable  22RV1  cell  lines  using  the  lentiviral  derived  tetracycline  inducible 
shRNA  knock  down  system.  We  were  unable  to  establish  an  inducible  SLC02B1  deficient  stable 
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cell  line  in  LNCaP,  thus,  we  used  a  transient-inducible  knock  down  of  SLC02B1  LNCaP  cells  in 


this  study  (shRNA-SLC02Bl).  After  successfully  knocking  down  SLC02B1  (Fig.  1A),  DHEAS 
uptake  was  substantially  decreased  to  -50%  in  22RV1  and  -70%  in  LNCaP  of  that  observed  in 
control  cells  (cells  transfected  with  scrambled  shRNA)  (Fig.  1C).  These  results  indicate  that 
SLC02B1  plays  an  essential  role  in  DHEAS  import  into  PC  cells  (Fig.  1C).  More  importantly, 
knocking  down  SLC02B1  abolished  atorvastatin’s  inhibition  of  DHEAS  uptake;  further 
suggesting  that  atorvastatin  competes  with  DHEAS  for  binding  to  their  transporter,  SLC02B 1 . 

Inhibition  of  SLC02B1  mediated  DHEAS  uptake  by  atorvastatin  decreases  PC  cell 
proliferation 

To  further  support  our  hypothesis,  we  wanted  to  demonstrate  that  inhibition  of  cell  growth  was 
SLC02B 1  and  DHEAS  dependent.  Thus,  we  examined  the  impact  of  atorvastatin,  the  most 
commonly  used  statin  clinically,  and  DHEAS  on  tumor  proliferation  before  and  after  SLC02B 1 
was  knocked  down  using  the  WST-1  assay. 

2.5  pM  (or  a  higher  concentration)  atorvastatin  is  known  to  inhibit  LNCaP  cell  proliferation  and 
induces  autophagy.13  Since  atorvastatin  concentrations  in  patient  serum  range  from  5-270nM,14 
we  chose  200  nM  of  atorvastatin  for  our  cell  proliferation  assay.  80  nM  DHEAS  significantly 
increased  cell  proliferation  in  LNCaP  and  22RV 1  lines  that  were  maintained  in  androgen- 
depleted  medium.  Cell  numbers  nearly  doubled  by  day  6  for  LNCaP  (Fig.  2A).  At  day  6, 
DHEAS  induced  a  -  6-fold  increase  in  LNCaP  cell  number,  compared  to  an  only  -3  fold 
increase  in  the  absence  of  DHEAS  (p=  0.0003).  However,  treatment  with  atorvastatin  (200  nM) 
inhibited  this  DHEAS  induced  cell  proliferation.  Consistent  with  this  finding,  knocking  down 
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SLC02B1  abolished  DHEAS  induced  cell  proliferation  in  LNCaP  and  22RV1  (Fig.  2B). 
Furthermore,  treatment  with  atorvastatin  did  not  significantly  inhibit  the  growth  of  SFC02B1 
knocked-down  cells.  Taken  together,  our  data  demonstrate  that  atorvastatin  can  efficiently  block 
SFC02Bl-mediated  DHEAS  uptake  and  DHEAS  induced  cell  growth  in  androgen  dependent 
PC  cell  lines.  This  mechanism  may  explain  the  observation  that  statin  use  may  be  associated 
improved  clinical  outcomes  in  PC  and  led  us  to  pursue  the  current  clinical  study. 

Clinical  Cohort  Study 

Of  the  926  patients  included,  283  patients  (30.6%)  were  taking  a  statin  at  ADT  initiation.  The 
majority  (93%)  of  patients  remained  on  a  statin  at  the  time  of  progression  or  at  last  follow-up. 
Statin  use  generally  increased  over  time  from  1996-2013  (eFig  2).  Of  the  643  non-users,  72 
started  a  statin  while  on  ADT;  in  these  patients,  median  time  from  ADT  initiation  to  initial  statin 
use  was  24  months  (IQR  12  to  ~48  months). 

Patient  and  disease  characteristics  at  diagnosis  and  at  ADT  initiation  are  detailed  in  Table  1. 
Statin  users  tended  to  have  a  lower  median  PSA  both  at  diagnosis  (9.1  vs.  11.8)  and  at  ADT 
initiation  (10.3  vs.  12.5).  Duration  from  diagnosis  to  ADT  initiation  was  longer  in  statin  patients 
(3.9  vs.  2.3  years).  Users  were  more  likely  to  have  lower  stage  disease  (56%  vs.  44%  T1  disease) 
and  less  likely  to  have  de  novo  metastases  (11%  vs.  18%)  or  nodal  involvement  (5%  vs.  10%)  at 
diagnosis.  Statin  users  were  more  likely  to  have  undergone  local  therapy  or  to  have  received 
ADT  as  part  of  local  therapy  (33%  vs.  26%).  Patients  on  statins  were  less  likely  to  have 
metastases  at  ADT  initiation  (53%  versus  63%)  (p  values  <0.05,  Table  1). 
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At  the  time  of  data  capture,  70%  (n=644)  of  patients  had  progressed  on  ADT  by  PSA.  Median 


follow-up  was  5.8  years  (range:  0.1-15.9).  Median  TTP  on  ADT  for  all  patients  irrespective  of 
statin  use  was  20.3  months  (95%  Cl:  17.5,  23.6).  Statin  users  at  ADT  initiation  had  a 
significantly  longer  median  TTP  on  ADT  (27.5  vs.  17.4  months,  p=0.0005,  Fig  3A).  Since  some 
but  not  all  of  the  baseline  characteristics  favored  more  indolent  disease  in  the  statin  users,  we 
adjusted  for  these  potential  imbalances.  Importantly,  the  association  remained  statistically 
significant  after  adjusting  for  pre-defined  prognostic  clinical  factors  including  biopsy  Gleason 
score,  type  of  primary  therapy,  use  of  prior  ADT  in  conjunction  with  localized  therapy, 
metastatic  status  and  PSA  at  initiation  of  ADT  (adjusted  HR=0.83,  95%  Cl:  0.69,0.99,  Table 
2).  When  stratified  by  year  of  ADT  initiation  using  5-year  increments,  the  association  between 
statin  use  and  TTP  on  ADT  remained  significant  in  the  multivariable  model  (adjusted  HR=0.83, 
95%  Cl:  0.69,1.00).  Moreover,  the  association  between  statin  use  and  TTP  was  observed 
regardless  of  whether  patients  had  radiographic  evidence  of  metastatic  disease  compared  to 
biochemical  relapse  only  at  ADT  initiation  (HR=0.79,  95%  CI:0.58,  1.07  for  M0;  HR=0.84,  95% 
0:0.67,1.06  for  Ml,  p  for  interaction=0.72).  (Fig.  3B,  3C) 


Discussion 

Progression  to  CRPC  depends  in  part  on  the  continued  reliance  of  the  tumor  on  AR  signaling  and 
residual  androgens.  In  our  preclinical  studies,  we  demonstrated  that  DHEAS  and  various  statins 
compete  for  binding  to  the  transporter  SLC02B 1  and  that  treatment  with  statins  competitively 
inhibit  DHEAS  uptake.  We  demonstrated  that  the  adrenal  androgen  DHEAS,  an  important 
precursor  to  DHT,  stimulates  PC  cell  proliferation  and  that  a  statin  drug,  atorvastatin,  can 
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diminish  DHEAS-stimulated  proliferation.  These  findings  drove  us  to  query  whether  statin  use 
might  influence  TTP  on  ADT  in  patients  with  hormone-sensitive  PC.  Given  the  vagaries  of 
defining  progression  in  PC,  in  our  clinical  study,  we  narrowly  defined  TTP  by  PSA  alone,  which 
allows  us  to  use  PSA  as  a  pharmacodynamic  endpoint  of  androgen  action. 

Statins  are  administered  widely  in  the  United  States  for  their  clinically  meaningful  lipid  lowering 
properties.  Most  epidemiological  studies  have  shown  significant  associations  between  statin  use 
and  decreased  incidence  of  advanced  prostate  cancer,  risk  of  recurrence  after  local  treatment, 
mortality,  and  PSA  levels  relative  to  non-users.5’6'15”17  Conversely,  most  cohort  and  case- 
controlled  studies  have  shown  no  association  between  statins  and  overall  prostate  cancer 
risk.  ’  ’  ’  To  our  knowledge,  none  have  reported  on  the  impact  of  statin  use  on  TTP  on  ADT. 


Given  the  biologic  heterogeneity  observed  in  prostate  cancer,  what  becomes  clear  is  that  more 
important  than  overall  incidence  is  the  impact  of  statin  use  on  the  development  of  lethal  prostate 
cancer.  A  recent  meta-analysis  collated  data  from  27  observational  studies  encompassing 
nearly  2  million  patients  and  assessed  the  association  of  statin  use  on  the  risk  of  developing  PC.5 
The  pooled  analysis  revealed  a  7%  reduction  in  the  risk  of  developing  any  PC  (RR  0.93,  95%  Cl 
0.87-0.99,  p=0.03).  Seven  of  the  studies  specifically  assessed  the  association  between  statin  use 
and  the  risk  of  clinically  significant  or  advanced  PC.  All  but  one  showed  a  relative  risk  reduction 
ranging  from  7-49%  with  statin  use  (RR:  0.51-0.93).  As  reviewed  by  Mucci  and  Stampher,  there 
have  been  at  least  5  additional  published  studies  since  2012  that  demonstrate  an  inverse 
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association  between  statin  use  and  prostate  cancer  mortality." 
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Our  study  evaluated  the  impact  of  statin  use  in  a  more  advanced  patient  population,  all  of  whom 
had  either  biochemical  or  metastatic  recurrence  after  local  therapy  or  de  novo  metastatic  disease 
for  which  they  were  started  on  ADT.  We  found  that  statin  use  at  the  time  of  ADT  initiation  was 
associated  with  a  significant  increase  in  TTP  on  ADT  even  after  adjusting  for  established 
prognostic  factors  such  as  Gleason  score  at  biopsy,  type  of  primary  therapy  be  it  radical 
prostatectomy  or  radiation,  use  of  ADT  with  the  primary  therapy,  or  presence  of  metastases  and 
PSA  at  ADT  initiation.1112 

Statins  inhibit  3-hydroxy-3-methylglutaryl-coenzyme  A  (HMG-CoA)  reductase,  an  enzyme 
critical  in  cholesterol  biosynthesis.”  Cholesterol  is  an  essential  component  in  the  synthesis  of 
steroid  hormones  such  as  the  androgens  that  drive  PC.  Other  purported  mechanisms  by  which 
statins  may  exert  their  antitumor  effects  include  inhibition  of  cell  proliferation,  inflammation, 
angiogenesis,  invasion,  and  metastasis  as  well  as  induction  of  apoptosis  and  autophagy.5,22"26  A 
reduced  pool  of  cholesterol  building  blocks  may  stymie  cancer  growth  as  rapidly  proliferating 
cells  like  tumor  cells  require  high  levels  of  available  cholesterol.”'  Further,  inhibition  of  HMG- 
CoA  results  in  reduction  of  mevalonate  and  downstream  famesylpyrophosphate  and 
geranylpyrophosphate,  which  are  important  in  activating  Ras  and  Rho.”  ’  Diminished  Ras  and 
Rho  activity  may  inhibit  cellular  proliferation  while  stimulating  apoptosis.26,27 

Multiple  preclinical  models  have  assessed  the  direct  effects  of  statins  on  prostate  cancer.  Zheng 
and  colleagues  evaluated  atorvastatin  and  celecoxib  in  xenograft  models.”  They  observed  a 
reduction  in  PC  cell  growth  even  after  castration  and  a  delay  in  progression  to  androgen 
independence.  In  vitro ,  Murtola  et  al  evaluated  the  effects  of  simvastatin  using  cell  lines  ranging 
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from  normal  prostate  epithelium  to  resected  primary  tumors  to  more  advanced  CRPC-like  cell 
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lines  such  as  LNCaP  and  VCaP.”  They  found  that  simvastatin  inhibited  early  stage  cell  linesbut 
not  the  more  advanced  CRPC-like  subtypes.  ’  While  not  directly  assessing  mechanism,  meta¬ 
analysis  of  multiple  randomized  controlled  trials  and  individual  cohort  studies  in  humans  have 

9Q  TO 

shown  that  statins  lower  testosterone  levels.”  ’ 

Using  PC  cell  lines,  we  found  that  statins  and  DHEAS  compete  for  the  same  transporter, 
SLC02B1,  and  that  administration  of  statins  competitively  reduces  uptake  of  DHEAS  and 
subsequent  tumor  cell  proliferation  (atorvastatin).  Thus,  statins  effectively  decrease  the  available 
intratumoral  androgen  pool,  which  may  enhance  cancer  control.  These  cell  line  findings  provide 
a  plausible  mechanism  to  support  our  clinical  observations  of  increased  TTP  on  ADT  in  statin 
users. 

An  important  limitation  of  our  study  is  the  retrospective  nature  in  which  we  gathered  data  on 
statin  use.  As  such,  usage  information  was  missing  in  16%  of  the  patients  in  our  cohort.  Several 
clinical  factors  significantly  differed  in  statin  users  compared  to  non-users.  At  ADT  initiation, 
patients  on  statins  tended  to  have  less  advanced  PC  with  less  nodal  or  metastatic  disease.  In 
contrast,  statin  users  were  more  likely  to  have  received  prior  ADT  with  local  therapy,  which 
might  have  impacted  the  durability  of  their  subsequent  use.  Nonetheless,  we  adjusted  for  these 
imbalances  in  our  multivariate  analysis.  Another  confounding  variable,  which  we  did  not 
specifically  address,  was  the  contribution  of  different  statins  on  TTP.  It  is  conceivable  that 
differences  in  statin  potency  and/or  their  pharmacokinetic  properties”  could  have  influenced  the 
magnitude  of  the  results. 
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Conclusion 


We  have  shown  that  statins  compete  with  DHEAS  for  influx  by  SLC02B1,  which  may  decrease 
the  tumor’s  available  androgen  pool.  Clinically,  this  may  translate  to  improved  cancer  outcomes, 
which  we  observed  in  our  institutional  cohort  of  men  receiving  primary  ADT.  Even  when 
controlled  for  known  prognostic  factors,  men  on  statins  had  a  significantly  longer  TTP  on  ADT 
than  non-users.  The  mechanism  through  which  statins  exert  their  activity  in  PC  is  likely 
multifactorial  including  antiproliferative  and  proapoptotic  effects  but  most  plausibly  is  the 
reduction  in  the  tumor’s  androgen  stores  through  a  combination  of  decreased  availability  of  the 
cholesterol  precursor  required  for  cle  novo  synthesis  and  decreased  transport  of  existing  precursor 
androgens  like  DHEAS  via  competitive  binding  of  SLC02B1.  The  widespread  use  of  statins  and 
their  established  safety  profile  make  them  attractive  potential  anticancer  therapeutics  as 
adjuvants  to  cytotoxic  or  androgen  ablating  therapies  or  as  preventative  agents.  Ultimately,  these 
results  require  prospective  validation.  Over  10  prospective  trials  are  ongoing  or  maturing  that 
will  further  characterize  the  role  of  statins  as  anticancer  therapies  in  PC.  (clinicaltrials.gov 
accessed  11/22/14). 
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Figure  Legend 


Figure  1.  Statins  inhibit  DHEAS  uptake  by  SLC02B1  in  PC  cells.  A.  Relative  mRNA  levels 
in  PC  cells  before  and  after  SLC02B1  is  knocked  down.  B.  DHEAS  uptake  in  PC  cells  with  2.5 
pM  DHEAS  and  different  concentrations  of  statins  when  incubated  for  60  minutes.  Statistical 
analysis  was  performed  by  comparing  each  condition  with  the  DHEAS  2.5  pM  and  no  statin 
state  (NS)  except  when  indicated.  C.  Depicts  DHEAS  uptake  in  PC  cells  before  (scrambled 
shRNA  (scr))  and  after  (shRNA  2B1)  SLC02B1  is  knocked  down  when  incubated  with  2.5  pM 
DHEAS  and  100  pM  atorvastatin  (ATO)  for  10  and  60  minutes.  Statistical  analysis  was 
performed  by  comparing  each  condition  with  scrambled  shRNA  (scr)  after  10  min  with  DHEAS 
except  when  indicated. 

Figure  2.  Atorvastatin  (ATO)  decreases  DHEAS  induced  PC  cell  proliferation  through 
SLC02B1. 

Atorvastatin  (ATO)  decreases  DHEAS  induced  PC  cell  proliferation  by  competing  for  influx 
through  SLC02B1.  A.  WST-1  assay  of  cell  proliferation  in  PC  cell  lines.  B.  WST-1  assay  of  cell 
proliferation  in  PC  cell  lines  after  SLC02Blis  knocked  down.  For  proliferation  assays,  cells 
were  maintained  in  androgen  depleted  medium  with  doxycycline  for  2  days  followed  by  the 
addition  of  80  nM  DHEAS  and  200  nM  ATO  to  the  culture  medium  for  the  indicated  time. 
Relative  cell  numbers  were  calculated  as  percentages  of  the  cell  numbers  at  day  0  (100%). 


Figure  3:  Kaplan  Meier  analysis  of  TTP  on  ADT  according  to  statin  use  in  all  patients. 

Kaplan  Meier  plot  of  TTP  on  ADT  according  to  statin  use  (A)  In  all  patients  (B)  In  patients 
with  no  visible  radiographic  metastasis  but  biochemical  (PSA)  failure  (M0)  and  (C)  in  patients 
with  metastasis  (Ml)  at  ADT  initiation.  (A)  Across  the  entire  cohort,  patients  on  statins  at  the 
time  of  ADT  initiation  had  a  significantly  longer  median  TTP  on  ADT  (27.5  vs.  17.4  months, 
p=0.0005).  The  association  between  statin  use  and  TTP  was  observed  regardless  of  whether 
patients  had  evidence  of  metastatic  disease  (C)  compared  to  biochemical  relapse  (B)  only  at 
ADT  initiation  (HR=0.79,  95%  CI:0.58,  1.07  for  M0;  HR=0.84,  95%  CI:0.67,1.06  for  Ml,  p  for 
interaction^  .72). 
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Supplementary  Figures: 

eFigure  1:  The  impact  of  statins  on  DHEAS  uptake  by  prostate  cancer  cells. 

DHEAS  uptake  in  PC  cells  with  100  pM  DHEAS  and  different  concentrations  of  statins  when 
incubated  for  60  minutes.  Statistical  analysis  was  performed  by  comparing  each  condition  with 
the  DHEAS  2.5  pM  and  no  statin  state  (NS). 

eFigure  2:  Statin  use  by  year  of  ADT  initiation.  Statin  use  tended  to  increase  over  time  with 
significantly  more  users  after  2006  compared  to  prior  to  2006. 
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Figure  1:  Statins  inhibit  DHEAS  uptake  by  SLC02B1  in  prostate  cancer  cells. 
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Figure  2:  Atorvastatin  (ATO)  decreases  DHEAS  induced  PC  cell  proliferation  through 
SLC02B1. 


LNCaP 


NT 

DHEAS 

ATO 

ATO  +  DHEAS 


day 


22RV1 


day 


day 


39 


Figure  3:  Kaplan  Meier  analysis  of  TTP  on  ADT  according  to  statin  use  in  all  patients. 
Kaplan  Meier  plot  of  TTP  on  ADT  according  to  statin  use  (A)  In  all  patients  (B)  In  patients 
with  no  visible  radiographic  metastasis  but  biochemical  (PSA)  failure  (MO)  and  (C)  in 
patients  with  metastasis  (Ml)  at  ADT  initiation. 
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Table  1  Patient  and  disease  characteristics  at  diagnosis  and  at  ADT  initiation  by  statin  use 
in  the  GELB  cohort  (N=926) 


N 

All 

% /Range 

Statin  Use 

No  (n=643) 

N  %  /Range 

Yes  (n=283 

N  % /Range 

p-value 

At  Diagnosis 

Age,  years  (median,  IQR) 

881 

61(55,67) 

616 

60(55,66) 

265 

62(56,67) 

0.018 

PSA,  ng/mL  (median,  IQR) 

813 

10.7(6,29) 

570 

11.8(6.3,40) 

243 

9.1(5.5,17) 

0.0003 

Clinical  T  stage 

T1 

442 

47.7 

283 

44.0 

159 

56.2 

0.005 

T2 

206 

22.2 

149 

23.2 

57 

20.1 

T3-4 

48 

5.2 

38 

5.9 

10 

3.5 

Tx/Unknown 

230 

24.9 

173 

26.9 

57 

20.1 

Clinical  M  stage 

MO 

338 

36.5 

235 

36.5 

103 

36.4 

0.013 

Ml 

148 

16.0 

117 

18.2 

31 

11.0 

Mx/Unknown 

440 

47.5 

291 

45.3 

149 

52.7 

Clinical  N  stage 

NO 

337 

36.4 

235 

36.5 

102 

36.0 

0.031 

N1 

78 

8.4 

64 

10.0 

14 

5.0 

Nx/Unknown 

511 

55.2 

344 

53.5 

167 

59.0 

Biopsy  Gleason  Score 

<6 

134 

14.5 

85 

13.2 

49 

17.3 

0.107 

7 

310 

33.5 

207 

32.2 

103 

36.4 
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Statin  Use 

All 

No  (n=643) 

Yes  (n=283 

N 

% /Range 

N 

% /Range 

N 

% /Range 

p-value 

At  Diagnosis 

Age,  years  (median,  IQR) 

881 

61(55,67) 

616 

60(55,66) 

265 

62(56,67) 

0.018 

PSA,  ng/mL  (median,  IQR) 

813 

10.7(6,29) 

570 

11.8(6.3,40) 

243 

9.1(5.5,17) 

0.0003 

>8 

376 

40.6 

272 

42.3 

104 

36.7 

Unknown 

106 

11.4 

79 

12.3 

27 

9.5 

Type  of  local  therapy 

RP  +/-RT 

388 

41.9 

268 

41.7 

120 

42.4 

<.0001 

RT  only/other 

285 

30.8 

172 

26.7 

113 

39.9 

None 

253 

27.3 

203 

31.6 

50 

17.7 

ADT  as  part  of  local  therapy 

No 

668 

72.1 

478 

74.3 

190 

67.1 

0.024 

Yes 

258 

27.9 

165 

25.7 

93 

32.9 

At  ADT  Initiation 

PSA,  ng/mL  (median,  IQR) 

849 

11.8(4.4,45) 

580 

12.5(4.4,59.1 

) 

269 

10.3(4.6,28.2 

) 

0.038 

Years  from  diagnosis  to  ADT 

881 

2.88(0.16,6.1 

616 

2.33(0.1,5.6) 

265 

3.85(1.2,7.8) 

<.0001 

initiation  (median,  IQR) 

) 

Metastases 

No 

372 

40.2 

239 

37.2 

133 

47.0 

0.005 

Yes 

554 

59.8 

404 

62.8 

150 

53.0 
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Statin  Use 

All 

No  (n=643) 

Yes  (n=283 

N 

% /Range 

N 

% /Range 

N 

% /Range 

p-value 

At  Diagnosis 

Age,  years  (median,  IQR) 

881 

61(55,67) 

616 

60(55,66) 

265 

62(56,67) 

0.018 

PSA,  ng/mL  (median,  IQR) 

813 

10.7(6,29) 

570 

11.8(6.3,40) 

243 

9.1(5.5,17) 

0.0003 

Concomitant  use,  prior  to  progression, 
of 

5 -alpha  reductase  inhibitor 

32 

3.5 

25 

3.9 

7 

2.5 

0.278 

Antiandrogen 

658 

71.1 

459 

71.4 

199 

70.3 

0.742 

Prior  chemotherapy 

72 

7.8 

40 

6.2 

32 

11.3 

0.008 
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Table  2:  Association  of  statin  use  with  TTP  on  ADT 


Statin  users 
(N=283) 

Non  users 
(N=643) 

P  value 

N  (%)  events 

164(58) 

480(75) 

Median  TTP  (95%  Cl), 
months 

27.5  (21.1,37.7) 

17.4(14.9,21.1) 

0.0005 

Univariable  HR  (95%  Cl) 

0.73  (0.61,0.87) 

Reference 

0.0005 

Multivariable  1  HR* 

0.83  (0.69,0.99) 

Reference 

0.0392 

Multivariable  2  HR  ** 

0.83  (0.69,1.00) 

Reference 

*The  multivariable  model  1  was  adjusted  for  biopsy  Gleason  score  (<7,  7,  >7),  type  of  primary 


therapy,  use  of  prior  ADT  in  conjunction  with  local  therapy,  metastatic  status  and  PSA  at 
initiation  of  ADT  (<  10, 10-20, >20ng/mL).  An  unknown  group  was  included  if  there  were 
missing  values  for  a  factor. 

**The  multivariable  model  2  was  a  stratified  Cox  regression,  adjusting  for  the  same  variables  as 
in  model  1  and  also  stratified  by  year  of  ADT  initiation  using  5-year  increments. 
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Supplementary  Data 


Supplementary  Methods: 
shRNA 

We  constructed  the  inducible  shSLC02Bl  RNA  expressing  plasmid  using  the  “all-in-one” 
pLKO-Tet-On  lentiviral  vector.  The  sequence  of  SLC02B1  shRNA  oligos:  forward  5'- 
ccggGAGGAGAGGGTTTGCTAATCTctcgagAGATTAGCAAACCCTCTCCTCttttt-3'  and 
reverse  5'- 

aattaaaaaGAGGAGAGGGTTTGCTAATCTctcgagAGATTAGCAAACCCTCTCCTC-3'. 

Quantitative  RT-PCR 

The  following  primers  were  used  to  measure  the  targeted  gene  expression  level:  SLC02B1, 
forward  5'-GTTTCGGCGAAAGGTCTTAGCAG-3',  reverse  5'- 
CCATCCTGCTTCTTCGTGGACT-3'  and  GAPDH,  forward  5'- 

CAGCCTCAAGATCATCAGCA-3',  reverse  5'-GTCTTCTGGGTGGCAGTGAT-3'.  The  primer 
pair  efficiencies  for  SLC02B1  and  GAPDH  were  calculated  according  to  the  linear  regression  of 
the  log  serial  (8  x)  dilution  factor  versus  the  correspondent  threshold  cycle  (Ct)  values.  The  slope 
of  the  regression  line  was  used  to  calculate  the  primer  pair  efficiencies  (E),  i.e.,  E  =  IOM/slope|. 

The  primer  pairs  with  efficiencies  close  to  2  were  used. 
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DHEAS  (pm ol /mg  protein) 


Supplementary  Figure  1:  The  impact  of  statins  on  DHEAS  uptake  by  prostate  cancer  cells. 
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Supplementary  Figure  2:  Statin  use  by  year  of  ADT  initiation 
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Purpose 

We  sought  to  validate  the  association  of  three  previously  demonstrated  germline  variants  in 
SLC02B1  with  time  to  progression  (TTP)  on  androgen  deprivation  therapy  (ADT)  and  to 
evaluate  if  these  genetic  variants  were  associated  with  overall  survival  (OS)  for  prostate  cancer 
(PC). 

Patients  and  Methods 

The  three  single  nucleotide  polymorphisms  (SNPs)  were  genotyped  in  an  independent  validation 
cohort  of  616  patients  with  PC  treated  with  ADT  at  our  institution  from  1996  to  2013. 
Multivariable  Cox  regression  adjusting  for  known  prognostic  factors  estimated  the  association  of 
genetic  variants  with  TTP  on  ADT  and  OS.  The  expression  of  SLC02B 1  was  examined  in 
prostatectomy  samples.  The  impact  of  SLC02B1  expression  level  on  uptake  of  DHEAS  was 
evaluated  in  cell  lines. 

Results 

The  association  between  the  exonic  SNP  rsl2422149  and  TTP  on  ADT  was  confirmed  in 
univariable  fP=.019)  and  multivariable  (adjusted  HR=  1.31,  95%  Cl:  (1.00,  1.72)  for  GG  vs. 
AA/AG,  P=.049)  analyses.  Since  OS  had  not  been  previously  evaluated,  we  examined  the 
association  in  the  combined  initial  and  validation  cohorts  (n=1094).  The  intronic  SNP  rs  1077858 
was  associated  with  OS  in  both  univariable  (P=.009,  Bonferroni  adjusted  P=. 027)  and 
multivariable  (adjusted  HR=  1.35,  95%  Cl:  (1.07,  1.71)  for  GG  vs.  AA/AG,  P=.012)  analyses. 
SLC02B1  expression  in  normal  prostate  tissue  and  in  22RV1  cells  carrying  the  major  allele  of 
SNP  rs  1077858  was  significantly  lower  than  those  carrying  the  risk  allele.  In  vitro,  we  showed 
SLC02B1  expression  levels  correlated  with  DHEAS  uptake  by  PC  cells. 

Conclusion 
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The  association  of  the  SNP  rs  1077858  with  OS  may  be  in  part  due  to  differential  SLC02B1 
expression  and  consequent  higher  uptake  of  DHEAS  and  subsequent  resistance  to  ADT,  which  in 
turn  might  contribute  to  shortened  survival. 
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INTRODUCTION 


Androgen  deprivation  therapy  (ADT)  is  the  primary  treatment  for  newly  diagnosed 
metastatic  prostate  cancer  (PC)  by  suppressing  testicular  androgen  production  and  diminishing 
AR  activation.  ’  Most  patients  will  eventually  progress  to  castration  resistant  PC  (CRPC)  which 
is  generally  lethal.3  6  The  variability  of  time  to  castration-resistant  disease  may  be  partly  due  to 
inherited  factors  that  result  in  increased  uptake  or  utilization  of  available  circulating  androgens. 

Despite  treatment  with  ADT,  persistent  AR  activation  stimulated  by  residual  androgen 
remains  an  integral  factor  that  drives  progression  to  castration  resistance.  '  These  residual 
androgens  are  derived  either  from  increased  de  novo  androgen  synthesis  or  by  uptake  and 
conversion  of  circulating  adrenal  androgens.6'9  One  such  adrenal  androgen, 
dehydroepiandrosterone  (DHEA)  and  its  sulfated  form,  DHEAS  are  secreted  in  large  amounts  by 
the  adrenal  cortex.  In  the  tumor  microenvironment,  DHEAS  is  converted  to  testosterone,  and 
subsequently  to  dihydrotestosterone  (DHT),  which  is  the  primary  androgen  utilized  by  PC 
cells.10"12  Even  after  ADT,  in  humans,  DHEAS  persists  at  a  high  concentration  in  the  blood  ata 
range  of  1-20  pM. 13-15 

How  patients  utilize  available  DHEAS  may  impact  their  response  to  ADT.  The  super  family 
of  organic  anion  transporting  polypeptides  (OATPs)  encoded  by  the  Solute  Carrier  Organ  Anion 
Transporter  Family  ( SLCO )  genes,  mediates  the  sodium-independent  uptake  of  a  wide  variety  of 
endogenous  compounds  and  drugs  into  cells.1617  SLC02B1,  one  of  the  SLCO  gene  family 
members,  mediates  the  transport  of  endogenous  steroid  conjugates,  such  as  DHEAS,  and  is 
expressed  in  a  broad  range  of  tissues  including  the  prostate.  CRPC  tumors  exhibit  higher 
SLC02B 1  expression  than  localized  PC  primary  tumors.  Epidemiologic  studies  have 
demonstrated  that  sequence  variations  of  SLC02B1  can  have  a  significant  impact  on  the  time  to 
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progression  (TTP)  to  CRPC  or  on  PC  survival."  '  In  our  previous  study,  we  genotyped  18 
single  nucleotide  polymorphisms  (SNPs)  of  SLC02B1  in  a  cohort  of  538  patients  with  PC  treated 
with  ADT.22  Three  SNPs  were  associated  with  TTP  on  ADT  (p<0.05):  an  exonic  SNP 
rsl2422149G>A  and  two  intronic  SNPs:  rsl789693A>T  and  rsl077858A>G.  The  differences  in 
median  TTP  for  each  of  those  three  identified  SNPs  were  approximately  10  months. 

In  this  report,  we  attempted  to  validate  whether  these  three  SLC02B1  SNPs  were  associated 
with  TTP  in  patients  on  ADT  in  an  independent  validation  cohort.  In  addition,  we  examined  the 
association  of  SLC02B1  SNPs  with  TTP  and  overall  survival  (OS)  in  the  combined  initial  and 
validation  cohorts  as  well  as  assessed  the  impact  of  metastatic  disease  status  at  ADT  initiation. 


PATIENTS  AND  METHODS 

All  cohorts  were  identified  from  our  established,  institutional  Prostate  Clinical  Research 
Information  System  (CRIS)  database."  Patients  with  biochemical  recurrence  (MO)  or 
radiologic  ally  evident  metastatic  disease  (Ml)  who  had  received  ADT  for  hormone  sensitive 
prostate  cancer  were  identified.  We  started  to  register  patients  to  CRIS  in  2001  and  all  patients 
were  consented  to  an  IRB-approved  protocol  that  collects  clinical  and  genomic  data.  For  the 
validation  part,  we  used  a  cohort  of  616  patients  who  had  not  been  previously  analyzed 
(validation  cohort).  For  other  analyses,  we  used  the  combined  cohorts  (N=1094),  which  also 
include  478  patients  from  the  initial  cohort  .  The  clinical  data  collection  as  well  as  the  detailed 
patient  characteristics  are  described  in  the  Supplemental  Data  and  have  been  previously 
described.22 
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Statistical  methods  were  described  in  the  Supplemental  Data.  Details  of  SNP  genotyping, 
plasmid  preparation,  cell  cultures,  transfections,  DHEAS  uptake  assays  were  described  in  the 
previous  study.22 

RESULTS 

Patient  Characteristics 

The  clinicodemographic,  disease,  and  past  treatment  characteristics  of  the  initial  (N=478)  and 
validation  ADT  cohort  (n=616)  are  presented  in  Supplementary  Table  1.  The  two  cohorts  had 
similar  disease  characteristics  at  diagnosis  and  at  ADT  initiation.  In  both  cohorts,  approximately 
70%  of  patients  had  received  a  local  therapy  (radical  prostatectomy  or  radiation  therapy)  and 
nearly  60%  had  metastases  at  the  time  of  ADT  initiation.  The  validation  cohort  had  a  lower  PSA 
at  diagnosis  (median  9.9  vs.  14  ng/ml),  probably  reflecting  the  more  frequent  use  of  PSA 
screening  leading  to  an  earlier  diagnosis  of  PC  in  those  patients.  The  validation  cohort  more 
frequently  received  intermittent  ADT  (35.9%  vs  18.8%). 

Validation  of  the  Association  of  SLC02B1  SNPs  with  TTP  during  ADT  in  the  validation 
cohort  (N=616) 

At  the  time  of  data  retrieval,  66%  (n=408)  of  patients  in  the  validation  cohort  had  progressed 
on  ADT.  The  median  TTP  on  ADT  was  20.9  months  (95%  Cl:  18.0,  24.0)  and  the  median 
follow-up  time  was  4.2  years  (range:  0.1,  16.3). 

The  exonic  SNP,  rsl2422149,  in  SLC02B1  was  associated  with  TTP  on  ADT  in  univariable 
analysis  (median  TTP:  27.2  (95%  Cl:  18.9,48.9)  months  for  AA/AG  and  20.0  (16.5,23.0)  months 
for  GG,  P  =  .019)(Table  1  and  Supplementary  Fig  1A).  In  multivariable  analyses  adjusting  for 
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clinical  factors,  the  association  between  the  SNP  and  TTP  remained  significant  (HR=1.31,  95% 
Cl:  1.00,  1.72,  P=.049).  Thus,  the  association  of  rsl2422149  with  TTP  on  ADT  observed  in  the 
initial  cohort  was  validated  such  that  PC  patients  carrying  the  major  GG  genotype  exhibited  a 
shorter  TTP  on  ADT.  While  there  was  a  trend  toward  an  association  between  rs  1077858  with 
TTP  on  ADT,  the  two  intronic  SNPs  were  not  validated  (Table  1  and  Supplementary  Fig  IB,  1C) 

Association  of  SLC02B1  SNPs  with  TTP  during  ADT  in  the  Combined  ADT  Cohorts  Based 
on  Presence  of  Metastases  at  ADT  Initiation  and  prior  ADT  use. 

Radiographically  evident  metastatic  PC  (Ml)  carries  a  worse  prognosis  than  biochemically 
recurrent  disease  (PSA  only  or  MO).  Thus,  we  explored  the  correlation  between  the  SLC02B1 
SNPs  and  TTP  in  the  combined  initial  and  validation  patients  and  stratified  by  metastatic  disease 
status  at  time  of  ADT  initiation,  both  of  which  had  not  been  previously  evaluated.  In  the 
combined  cohort  (n=1094),  74%  of  patients  (n=811)  had  progressed  on  ADT.  Median  TTP  on 
ADT  was  18.9  months  (95%  Cl:  16.5,  21.1). 

For  the  exonic  rs  12422 149,  we  again  found  a  significant  association  with  TTP  on  ADT  in  the 
combined  cohort  (adjusted  HR=1.33,  95%CI:  1.10-1.60,  P-  .003).  When  stratified  by  metastatic 
disease  status,  the  association  remained  in  the  Ml  population  (adjusted  HR=1.65,  95%:  1.28- 
2.12)  but  was  not  seen  in  the  M0  population  (adjusted  HR=0.96,  95%  Cl:  0.72-1.27)  ( P 
interaction=  .006,  Table  2  and  Supplementary  Fig  2).  The  association  between  the  intronic  SNP 
rsl077858  and  TTP  was  of  borderline  significance  (P=  .075;  P=  .062  if  AA  and  AG  were 
combined).  No  association  was  found  for  rs  1789693  and  TTP  on  ADT  in  the  combined  cohort. 
Results  were  similarly  negative  in  the  M0  and  Ml  populations  (P  interaction  >0.5,  Table  2  and 
Supplementary  Fig  2). 
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The  associations  with  TTP  on  ADT  were  similar  among  patients  with  and  without  prior 
hormone  treatments  for  the  exonic  SNP  rsl2422149.  However,  for  the  intronic  SNP  rsl077858,  a 
significant  association  was  observed  only  in  patients  without  prior  hormone  treatment  (adjusted 
HR=1.32  (1.06,  1.65)),  but  not  in  patients  with  prior  hormone  treatment  (adjusted  HR=0.87  (0.55, 
1.37))  (Supplementary  Table  2). 

Association  of  SLC02B1  SNPs  with  OS  in  the  Combined  ADT  Cohorts 

Correlations  with  OS  had  not  been  explored  in  our  prior  work  and  thus  we  evaluated  this 
endpoint  in  the  combined  cohort.  Nearly  half  (49%,  n=537/1094)  of  the  patients  had  died  at  the 
time  of  data  collection.  The  median  OS  from  ADT  initiation  was  6.5  years  (95%  Cl:  6.0,  7.0)  and 
the  median  follow-up  time  was  6.5  years  (range  0.1,  16.3  years). 

There  was  no  statistically  significant  association  with  OS  for  the  exonic  SNP  rsl2422149  in 
either  univariable  or  multivariable  analyses  (Table  3  and  Supplementary  Fig  3A).  In  univariable 
analysis,  there  was  no  association  between  rs  1789693  and  OS  from  ADT  initiation  ( P=  .184, 
Bonferroni  adjusted  P=0.552),  but  patients  carrying  the  minor  allele  (AT  or  TT)  had  longer  OS 
in  multivariable  analysis  (HR=0.81  and  0.78,  respectively,  P=  .044)  (Table  3  and  Supplementary 
Fig  3B).  Importantly,  we  found  that  patients  having  the  minor  genotype  GG  in  the  intronic  SNP 
rs  1077858  had  a  shorter  OS  from  ADT  initiation  in  both  univariable  ( P=  .009,  Bonferroni 
adjusted  P= 0.027)  and  multivariable  analysis  (adjusted  HR=1. 35(1.07,  1.71),  P=  .012)  (Table  3 
and  Supplementary  Fig  3C).  The  median  OS  was  decreased  from  6.7  (95%  Cl:  6.2,  7.2)  to  5.2 
(95%  Cl:  4.3,  6.8)  years  with  the  AA/AG  vs.  GG  rsl077858  genotypes  respectively. 

When  stratified  by  metastatic  disease  status  (MO  vs.  Ml),  patients  carrying  the  exonic  SNP 
rs  12422 149  GG  genotype  had  a  trend  for  shorter  OS  only  in  the  Ml  population  and  at  a 
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borderline  significance,  whereas  this  difference  was  not  observed  in  the  MO  patients  (adjusted 
HR  for  GG:  0.82  (95%CI:  0.58-1.14)  for  MO;  1.32(0.97-1.80)  for  Ml,  P  interaction=  .038).  The 
association  of  OS  with  rsl789693  and  rsl077858  were  similar  between  M0  and  Ml  population 
(P  interaction  >0.9,  Table  3  and  Supplementary  Fig  4).  The  minor  rsl077858  GG  genotype 
portended  shorter  OS  from  ADT  initiation  in  both  groups  (adjusted  HR=1.36  (95%CI:  0.93,  1.99) 
and  1.35  (95%CI:  1.00,  1.82)  respectively)  (Table  3  and  Supplementary  Fig  4). 

Impact  of  the  Intronic  SNP  rsl077858  on  Expression  of  SLC02B1 

In  our  previous  study,  we  showed  that  different  SLC02B1  variants  corresponding  to  the 
exonic  SNP  rsl2422149G>A  (Arg312Gln)  influenced  the  efficiency  of  DHEAS  transport  into 
PC  cells.22  Given  the  fact  that  the  intronic  SNP  rsl077858  was  associated  with  OS,  we  studied 
the  possible  mechanism  by  which  the  SNP  might  affect  OS  from  ADT  initiation.  Using  80 
normal  prostate  tissue  samples  available  in  our  tissue  repository  of  patients  known  to  have  PC, 
we  analyzed  SLC02B1  expression  (Fig  1A).  Results  demonstrated  that  the  different  rs  1077 85 8 
SNPs  were  associated  with  significantly  different  SLC02B1  expression  levels  (Ptrend=  .019)  with 
the  highest  SLC02B1  expression  in  patients  carrying  the  minor  GG  genotype.  To  prove  that 
different  rsl077858  SNPs  impact  the  expression  of  SLC02B1,  we  reconstituted  the  rsl077858 
SNP  variants  using  CRISPR  in  22RV1  cells  to  convert  the  original  A  allele  to  the  G  allele 
(Supplementary  data).  The  association  between  different  rs  1077858  SNPs  and  SLC02B1 
expression  levels  was  significant  (/V-nd=  .0073)  and,  once  again,  the  GG  genotype  in  22RV1 
cells  corresponded  to  the  highest  SLC02B1  expression  level  (Fig  IB).  These  results  suggest  that 
the  association  of  rs  1077858  on  OS  was  possibly  through  changes  in  SLC02B1  expression 
levels  in  these  patients.  As  a  result,  alleles  corresponding  to  higher  SLC02B1  expression  may 
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allow  greater  DHEAS  uptake  and  thereby  induce  shorter  OS  in  patients.  To  evaluate  this 


hypothesis,  we  examined  DHEAS  uptake  in  a  LNCaP  PC  cell  line  by  altering  the  expression 
level  of  SLC02B1  to  mimic  different  rsl077858  allele  conditions.  After  successfully  knocking 
down  or  overexpressing  SLC02B1  (Fig  1C),  we  found  that  DHEAS  uptake  activity  was 
dependent  on  the  expression  level  of  SLC02B1  and  that  higher  expression  of  SLC02B1 
transported  more  DHEAS  into  the  cells  (Fig  ID). 

DISCUSSION 

There  is  great  variability  in  the  durability  of  response  to  ADT,  and  multiple  mechanisms  are 
likely  operative.  Persistent  androgens,  especially  adrenal  androgen  precursors  such  as  DHEAS, 
play  an  important  role  in  the  evolution  of  CRPC  by  continuing  to  drive  the  AR  signaling 
pathway. 15,25,26  SLC02B1  is  an  active  androgen  transporter  that  pumps  androgens  and  androgen 
precursors  into  prostate  cells,  potentially  promoting  resistance  to  ADT7  '  ’  We  have  found 
that  SLC02B1  SNPs  are  significantly  associated  with  either  the  efficiency  of  DHEAS  transport 
or  the  expression  level  of  SLC02B 1  on  prostate  cells,  which  may  in  part  mechanistically 
determine  the  variability  in  outcomes  on  ADT. 

Our  initial  study  revealed  that  three  SLC02B1  SNPs  significantly  influenced  TTP  on  ADT  in 
PC  patients  with  approximately  10  months  differences  in  median  TTP.  “  In  addition,  an 
independent  study  of  the  SLC02B1  SNP  in  a  Japanese  population  also  found  that  rsl2422149 
associate  with  TTP  during  ADT.“  In  the  current  study,  we  attempted  to  further  validate  those 
findings  in  an  independent  cohort  of  patients  and  to  better  elucidate  the  mechanism  behind  these 
effects.  Only  the  exonic  SNP  rsl2422149  was  validated  for  the  association  with  TTP  on  ADT  in 
this  study,  such  that  PC  patients  who  carry  the  major  GG  genotype  for  rs  12422 149  may  derive 
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less  benefit  from  ADT.  Increased  uptake  in  the  poorer  risk  allele  carriers  may  be  associated  with 
greater  intracellular  androgen  stores  and  thus,  is  a  plausible  mechanism  of  resistance  to  ADT.  In 
addition,  in  subset  analysis,  the  association  of  rsl2422149  with  TTP  on  ADT  was  seen 
predominantly  in  Ml  patients,  further  suggesting  additional  factors  yet  to  be  elucidated  may 
participate  in  resistance  in  Ml  patients  compared  to  the  earlier  MO  disease  state. 

The  difference  in  median  TTP  between  the  AA  and  GG  rs  1077858  genotypes  in  the 
validation  cohort  was  6.8  months  (AA:  20.9  months  (95%  Cl:  17.0,  25.1)  vs.  GG:  14.1  months 
(95%  Cl:  9.7,  23.7)),  which  has  a  similar  trend  to  that  observed  in  our  initial  study.  This 
difference  did  not  persist  in  the  multivariable  model  after  adjusting  for  other  clinical  factors. 
False  discovery  is  certainly  a  possible  reason  for  our  inability  to  validate  an  association  between 
the  two  intronic  SNPs,  rsl789693  and  rsl077858,  and  TTP  on  ADT  in  the  current  study. 

To  our  knowledge,  the  current  study  is  the  first  one  to  demonstrate  that  an  inherited  variation 
in  the  SLC02B1  gene  was  significantly  associated  with  OS  after  ADT.  While  the  intronic  SNP 
rs  1077858  correlated  with  differences  in  TTP  in  the  initial  cohort,  the  current  study  should  be 
viewed  as  hypothesis  generating.  In  the  current  study,  it  was  found  to  be  associated  with  OS  in 
the  larger  combined  cohort  of  patients.  Importantly,  the  phenotype  of  this  functional  SNP  is 
independent  of  metastatic  disease  status  at  ADT  initiation,  indicating  that  the  decreased  OS 
cannot  be  attributed  to  lead  time  bias,  i.e.,  starting  ADT  later  when  patients  have 
radiographically  evident  metastases  compared  to  biochemical  recurrence  only.  Taking  it  a  step 
further  mechanistically,  we  demonstrated  that  genetic  variations  in  rsl077858  correlate  with 
variable  expression  levels  of  SLC02B 1  in  patient  prostate  tissue  samples  and  in  a  cell  line  based 
study.  As  such,  it  is  plausible  that  the  rsl077858  GG  variant  with  its  higher  level  of  SLC02B1 
expression  results  in  enhanced  DHEAS  uptake  and  increased  activation  of  the  AR  signaling 
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pathway  with  subsequent  shorter  OS  in  patients.  It  is  unclear  why  we  observed  an  OS  benefit  but 
no  definite  impact  on  TTP  with  rs  1077858  in  the  validation  cohort.  However,  there  was  a  median 
6  month  increase  in  TTP  which  was  consistent  with  an  improvement  in  outcomes  with  this  SNP 
but  this  observation  did  not  reach  statistical  significance  in  the  univariate  model  (TTP  on  ADT 
for  AA  and  GG  is  20.9  and  14.1  months,  respectively.  P=. 383). 

One  must  also  question  why  the  only  SNP  validated  to  influence  TTP  on  ADT,  the  exonic 
rsl2422149,  was  not  associated  with  OS.  It  is  worth  mentioning  that  although  not  associated 
with  OS,  GG  genotype  of  rsl2422149  had  shorter  TTP  and  OS  in  Ml  patients,  which  suggests 
that  rs  12422 149,  to  some  extent,  could  affect  the  TTP  and  OS  in  a  consistent  fashion.  Since  the 
clinical  evidence  for  a  significant  association  of  TTP  on  ADT  with  OS  has  not  been  established, 
and  the  exonic  SNP  rsl2422149  functionally  affecting  DHEAS  uptake  efficiency  has  a 
significant  association  with  TTP  on  ADT  but  not  with  OS  in  all  patients,  we  believe  that 
mechanisms  other  than  the  expression  level  of  SLC02B 1  may  be  involved  in  the  association  of 
rsl077858  with  OS.  For  example,  variability  in  DHEAS  uptake  may  be  an  important  driver  in 
earlier  stages  of  metastatic  disease  but  non-androgen  pathways  may  induce  resistance  in  later 
stages.  It  is  also  possible  that  unidentified  regulatory  or  functional  SNPs  crosstalk  with  the 
known  functional  SNPs  in  different  disease  stages  in  this  study.  For  example,  Heam  et.  al.  have 
found  that  men  inheriting  the  variant  HSD3B1  (1245C)  allele  that  enhances  DHT  synthesis 
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exhibit  resistance  to  ADT,  as  manifested  by  worse  clinical  outcomes  (shorter  survival  time).“ 
Therefore,  those  unidentified  factors,  such  as  SNP  in  HSD3B1  (A1245C),  need  to  be  considered 
in  future  analysis  to  further  characterize  the  association  of  SNPs  in  SLC02B1  with  TTP  and  OS. 
Finally  and  most  plausibly  with  respect  to  finding  a  TTP  but  not  OS  difference  is  heterogeneity 
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in  terms  of  access  and  administration  of  the  now  5  therapies  approved  for  CRPC  that  improve 
survival  including  enzalutamide,  abiraterone,  cabazitaxel,  radium-223,  and  sipuleucel-T.29'35 

Additional  potential  confounders  of  our  results  include  differences  between  the  initial  and 
validation  cohorts  in  terms  of  patient  and  disease  characteristics  at  diagnosis  and  at  ADT 
initiation.  Patients  from  the  validation  cohort  were  more  likely  to  have  a  lower  PSA  at  diagnosis 
and  to  have  received  hormonal  therapy  as  part  of  local  therapy  (Supplementary  Table  1).  Thus, 
the  validation  cohort  could  have  been  biased  toward  living  longer  in  one  respect  but  this  could  be 
countered  by  the  higher  use  of  prior  ADT,  which  would  portend  a  shorter  time  on  ADT.  When 
prior  use  of  ADT  was  added  to  the  multivariable  model,  it  was  not  found  to  impact  the  result. 
Finally,  differences  in  access  and  administration  of  subsequent  survival  improving  therapies  such 
as  enzalutamide  and  abiraterone  may  have  also  confounded  the  OS  results.  Meanwhile,  we  found 
that  the  exonic  SNP  rs  12422 149  effect  on  TTP  was  not  dependent  on  prior  ADT  but  the  effect  of 
the  intronic  SNP  rs  1077858  was  dependent  on  prior  ADT  (Supplementary  Table  2).  Thus,  future 
analyses  should  be  performed  stratifying  by  disease  stage,  use  of  ADT  with  local  therapy,  and 
subsequent  therapies  for  metastatic  disease  to  disclose  the  unknown  factors  that  affect  the 
different  SNP’s  functions  at  different  stages  of  the  disease  and  result  from  prior  ADT  use. 

In  conclusion,  germline  variants  within  SLC02B1  appear  to  modulate  function  or  expression 
of  SLC02B 1  in  vitro,  which  consequently  may  affect  the  uptake  of  DHEAS  and  impact  TTP  and 
OS  in  PC  patients.  These  findings  support  that  SLC02B1  variants  are  prognostic  markers  of 
durability  of  response  to  ADT.  Further,  SLC02B1  may  serve  as  a  target  for  novel  agents  and 
strategies.  Ultimately,  if  prospectively  validated,  genotyping  of  SLC02B1  for  patients  may  help 
tailor  more  effective  therapies  to  the  individual  patient. 
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Figure  Legends 

Fig  1.  Impact  of  the  intronic  SNP  rsl077858  on  the  expression  of  SLC02B 1  which  affects  the 
sulfated  dehydroepiandrosterone  (DHEAS)  uptake.  (A)  Relative  SLC02B1  mRNA  levels 
compared  to  GAPDH  in  normal  prostate  tissue  in  primary  prostate  cancer  patients.  Significant 
differences  were  observed  among  three  groups,  /Vend  <0.05.  Lines  in  the  figure  represent  the 
median  with  interquartile  range.  (B)  Relative  SLC02B1  mRNA  levels  compared  to  GAPDH  in 
22RV 1  cells  with  different  alleles  after  CRISPR  reconstitution.  Significant  difference  were 
observed  among  three  groups,  /VmKO.0 1 .  Lines  in  the  figure  represent  the  median  with 
interquartile  range.  (C)  SLC02B1  mRNA  levels  in  mock  transfected  LNCaP  cells  and  LNCaP 
cells  transfected  with  different  SLC02B1  variants.  Values  represent  the  fold  differences  relative 
to  those  in  mock  transfected  cells,  which  were  set  as  1.0.WT,  LNCaP  mock  transfected  with 
either  siRNA  control  or  pCMV6-XL4  vector;  siSLC02Bl,  LNCaP  transiently  transfected  with 
siRNA  targeting  SLC02B1  (-p  represents  the  SMARTpool  siRNA  from  GE  Healthcare;  -1,-2 
and  -3  represent  three  unique  27mer  siRNAs  from  Origene,  respectively.  See  Supplementary 
Data  for  details);  oeSLC02Bl,  LNCaP  transiently  transfected  with  pCMV6-XL4-SLC02Bl.  All 
experiments  were  repeated  in  triplex.  Statistical  analysis  (unpaired  t-test)  was  performed  by 
comparing  each  condition  with  the  WT  condition.  All  mRNA  levels  were  analyzed  by 
quantitative  reverse  transcriptase  polymerase  chain  reaction  and  normalized  by  the  expression 
level  of  glyceraldehyde  3-phosphate  dehydrogenase.  (D)  DHEAS  uptake  in  LNCaP  cells.  Cells 
were  incubated  with  2.5  pM  DHEAS  for  10  minutes  after  cells  have  been  transfected  with 
siRNA  or  pCMV6-XL4-SLC02Bl  for  3  days.  All  experiments  were  repeated  in  triplex. 
Statistical  analysis  (unpaired  t-test)  was  performed  by  comparing  each  condition  with  the  WT 
condition. 
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Table  1.  SLC02B1  Genotype  Distributions  and  Their  Association  with  TTP  on  ADT  in 

Validation  Cohort 


Univariate  Model 

Multivariate  Model** 

Genotype 

N(%) 

Median  TTP 
(95%  Cl), 
months 

Log-rank 

P  value 

HR 

(95%  Cl) 

Adjusted  HR 
(95%  Cl) 

P  value 

rs  12422 149 

AA/AG 

110(18) 

27.2(18.9,48.9) 

0.019 

1.00  (reference) 

1 .00  (reference) 

0.049 

GG 

506(82) 

20.0(16.5,23.0) 

1.39(1.05,1.80) 

1.31(1.00,1.72) 

rsl789693 

AA 

238(39) 

21.6(17.5,34.1) 

0.176 

1 .00  (reference) 

1.00  (reference) 

0.409 

AT 

281(46) 

19.4(15.0,23.4) 

1.22(0.99,1.51) 

1.15(0.93,1.43) 

TT 

90(15) 

25.2(18.1,38.1) 

1.07(0.79,1.45) 

1.02(0.75,1.39) 

rsl077858 

AA 

288(47) 

20.9(17.0,25.1) 

0.383 

1 .00  (reference) 

1.00  (reference) 

0.880 

AG 

242(40) 

23.0(16.6,30.0) 

(0.185)* 

0.96(0.77,1.18) 

0.98(0.79,1.22) 

(0.634)* 

GG 

81(13) 

14.1(9.7,23.7) 

1.19(0.88,1.59) 

1.06(0.79,1.44) 

*If  A  A  and  AG  were  combined. 

**  Adjusting  for  biopsy  Gleason  score,  type  of  primary  therapy,  use  of  prior  ADT  in 
conjunction  with  local  therapy,  metastatic  status  and  PSA  at  initiation  of  ADT. 
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Table  2.  Association  of  SLC02B1  Genotype  with  TTP  on  ADT  in  the  Combined  (Initial  plus 
Validation)  Cohort  and  Stratified  by  Metastatic  Disease  Status  at  ADT  Initiation 


*  Adjusting  for  biopsy  Gleason  score,  type  of  primary  therapy,  use  of  prior  ADT  in 
conjunction  with  local  therapy,  metastatic  status  and  PSA  at  initiation  of  ADT. 

**  From  multivariable  Cox  regression  models,  with  the  interaction  of  SNP  and  metastatic 
disease  and  simultaneously  adjusted  for  other  clinical  variables.  P  value  for  interaction  is  to 
test  whether  the  association  of  SNP  with  TTP  differs  by  metastatic  disease  status. 
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Table  3.  Association  of  SLC02B1  Genotype  with  OS  from  the  Initiation  of  ADT  in  the 
Combined  (Initial  plus  Validation)  Cohort  and  Stratified  by  Metastatic  Disease  Status  at 

ADT  Initiation 


Combined  Cohort 

By  Metastatic  Disease  Status*** 

Univariate  Model 

Multivariate 

M0 

Ml 

Model* ** 

* 

P 

Log- 

valu 

Genoty 

Median 

rank 

e 

pe 

N(%) 

OS 

P  value 

HR  (95% 
Cl) 

Adjusted 

P 

Adjusted 

Adjusted 

(inte 

(95% 

(Bonferr 

HR  (95% 

valu 

N 

HR  (95% 

N 

HR  (95% 

r- 

Cl), 

oni 

Cl) 

e 

Cl) 

Cl) 

actio 

years 

adjusted 

) 

n) 

rsl242 

2149 

AA/AG 

197(1 

7. 1(6.0, 

0.373 

1 .00(refere 

1.00(refere 

0.50 

93 

1.00(refere 

10 

1 .00(refere 

0.03 

8) 

8.1) 

(0.999) 

nee) 

nee) 

6 

nee) 

4 

nee) 

8 

GG 

894(8 

6. 4(5. 9, 

1.11(0.89,1 

1.08(0.86,1 

36 

0.82(0.58,1 

52 

1.32(0.97,1 

2) 

6.9) 

.39) 

.36) 

6 

.14) 

8 

.80) 

rsl789 

693 

AA 

427(3 

6. 1(5.2, 

0.184 

1 .00(refere 

1 .00(refere 

0.04 

18 

1.00(refere 

24 

1 .00(refere 

0.99 

9) 

6.7) 

(0.552) 

nee) 

nee) 

4 

6 

nee) 

1 

nee) 

9 

AT 

486(4 

7.2(6. 1, 

0.84(0.70,1 

0.81(0.67,0 

20 

0.81(0.60,1 

28 

0.81(0.64,1 

5) 

7.7) 

.01) 

.98) 

0 

.11) 

6 

.03) 

TT 

169(1 

6. 2(5.0, 

0.89(0.69,1 

0.78(0.60,1 

68 

0.78(0.50,1 

10 

0.77(0.57,1 

6) 

7.8) 

.14) 

.00) 

.21) 

1 

.06) 

rsl077 

858 

AA/AG 

939(8 

6. 7(6. 2, 

0.009 

1 .00(refere 

1 .00(refere 

0.01 

39 

1.00(refere 

54 

1 .00(refere 

0.97 

* 

6) 

7.2) 

(0.027) 

nee) 

nee) 

2 

7 

nee) 

2 

nee) 

0 

GG 

147(1 

5. 2(4.3, 

1.36(1.08,1 

1.35(1.07,1 

60 

1.36(0.93,1 

87 

1.35(1.00,1 

4) 

6.8) 

.71) 

.71) 

.99) 

.82) 

*Distributions  of  OS  were  similar  between  AA  and  AG  group,  thus  were  combined  (median 
6.8  and  6.5  years,  respectively) 

**Adjusting  for  biopsy  Gleason  score,  type  of  primary  therapy,  use  of  prior  ADT  in 
conjunction  with  local  therapy,  metastatic  status,  PSA  at  initiation  of  ADT  and  age  at  ADT 
initiation. 

***From  multivariable  Cox  regression  models,  with  the  interaction  of  SNP  and  metastatic 
disease  and  simultaneously  adjusted  for  other  clinical  variables.  P  value  for  interaction  is  to 
test  whether  the  association  of  SNP  with  TTP  differs  by  metastatic  disease  status. 
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SUPPLEMENTARY  DATA 


Patients  and  methods 

Briefly,  the  initial  ADT  cohort  (N=595)  was  identified  in  June  2006  and  the  patients’ 
treatments  and  outcome  data  were  updated  in  2012.  Of  these,  27  patients  were  excluded 
due  to  insufficient  follow-up  data  after  ADT  initiation  and  77  were  excluded  because 
they  started  ADT  prior  to  1996.  This  year  was  used  as  a  cutoff  because  patients  who  had 
lived  that  long  (i.e.,  at  least  6  yrs  from  time  of  ADT  initiation)  were  excellent  responders 
and  inclusion  would  have  created  the  potential  for  selection  bias.  Of  the  remaining  491 
patients,  478  were  genotyped  for  the  three  SLC02B1  SNPs.  In  the  validation  cohort,  758 
patients  were  identified  and  clinical  data  for  this  analysis  was  retrieved  from  the  CRIS 
database  in  2013.  Of  these,  104  patients  were  excluded  due  to  insufficient  follow-up  data 
after  ADT  initiation  and  1 1  were  excluded  because  they  started  ADT  prior  to  1996.  Of 
the  remaining  643  patients,  616  were  genotyped  for  the  three  SLC02B1  SNPs. 

TTP  was  defined  as  two  rises  in  prostate-specific  antigen  (PSA)  (at  least  1  week  apart) 
while  receiving  ADT.  The  first  rise  was  required  to  be  greater  than  the  nadir  PSA  value 
plus  0.02  ng/ml.  Initiation  of  a  second  treatment  for  rising  PSA  before  fulfillment  of  the 
definition  of  progression  was  also  considered  as  a  progression  event;  the  date  of  starting 
the  secondary  treatment  was  denoted  as  the  date  of  progression.  TTP  during  ADT  was 
defined  as  the  duration  of  time  from  ADT  initiation  to  the  date  of  ADT  progression  or 
the  date  of  initiation  of  secondary  therapy.  Among  patients  who  did  not  progress,  they 
were  censored  as  of  the  date  of  their  last  known  progression-free  visit  or  the  date  of  death 
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in  those  who  died  without  progression.  OS  was  defined  as  the  period  from  ADT 
initiation  to  patient  death  or  was  censored  at  date  of  last  follow  up. 


Statistical  Methods 

SNPs  were  analyzed  as  three  distinct  genotype  groups,  with  the  exception  of  the 
exonic  SNP  rs  12422 149  analysis  wherein  the  AA  and  AG  were  combined  because  only  1% 
of  the  population  carries  the  AA  genotype.  AA  and  AG  were  also  combined  for  the 
intronic  rs  1077858  SNP  during  the  analysis  of  their  association  with  OS  since  the 
distributions  of  OS  were  similar  between  the  A  A  and  AG  groups. 

Progression  was  defined  as  a  minimum  of  two  rises  in  prostate-specific  antigen  (PSA), 
with  the  date  of  first  rise  (nadir  +  >0.02  ng/ml)  as  the  progression  date.  TTP  during  ADT 
was  defined  as  the  duration  of  time  from  ADT  initiation  to  the  date  of  ADT  progression 
or  the  date  of  initiation  of  a  secondary  therapy  (for  rising  PSA  before  fulfillment  of  the 
definition  of  progression),  or  was  censored  at  the  date  of  last  follow-up  visit  or  PSA 
value  among  patients  who  did  not  progress.  OS  was  defined  as  the  period  from  ADT 
initiation  to  patient  death,  or  it  was  censored  on  the  day  of  last  follow  up. 

The  distributions  of  TTP  and  OS  were  estimated  using  the  Kaplan-Meier  method,  with 
95%  confidence  intervals  (CIs).  The  association  between  SNPs  and  TTP  was  evaluated 
using  the  log-rank  test  or  the  Wald  chi-square  test  and  by  multivariable  Cox  regression 
model  adjusted  for  known  prognostic  factors  including  biopsy  Gleason  score,  type  of 
primary  therapy,  use  of  prior  ADT  in  conjunction  with  local  therapy,  metastatic  status 
and  PSA  at  ADT  initiation.  OS  was  analyzed  similarly,  but  age  at  ADT  initiation  was 
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also  included  as  a  covariate  in  the  multivariable  models.  To  estimate  whether  the 
association  of  SNPs  with  TTP  during  ADT  and  OS  differed  by  metastatic  disease  status 
at  ADT  initiation,  multivariable  Cox  regression  models  were  performed  with  the 
interaction  term  of  SNP  and  metastatic  disease  status,  and  simultaneously  adjusted  for 
the  other  clinical  variables. 

For  TTP  analysis,  there  was  no  correction  for  multiple  comparisons  since  the 
validation  part  mainly  aimed  to  reproduce  the  strength  of  association  previously 
identified  in  our  initial  work  where  we  had  applied  appropriate  correction  for  multiple 
comparisons.  For  the  test  of  OS,  Bonferroni  procedure  was  used  to  correct  multiple 
comparisons. 

For  in  vitro  cell  line  studies,  data  was  represented  as  means  ±  standard  deviation  of  at 
least  3  biological  repeats.  Comparison  between  two  independent  groups  was  performed 
by  an  unpaired  2-tailed  t  test.  P  <  .05  (two  sided)  was  considered  statistically  significant 
for  all  analyses. 


siRNA 

The  efficiency  of  knocking  down  SLC02B1  expression  was  assayed  after  siRNA 
transfection  for  48  hours  by  RT-PCR  (Forward  oligo  is 
GTTTCGGCGAAAGGTCTTAGCAG  and  reverse  oligo  is 

CCATCCTGCTTCTTCGTGGACT,  which  were  purchased  from  Origene).  siSLC02Bl- 
p  was  purchased  from  GE  Healthcare  as  a  SMARTpool  with  4  target  sequences 
(CAUCCAUGGCUGCGGGCAU;  GCCACCAGAUUGCGGGCAU ; 
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UCUCGGAGCCAUACCGCUA;  AUAAUGACCUGCUCCGAAA).  siSLC02Bl-l,  -2, 
AND  -3  were  unique  25mer  siRNA  duplexes  and  purchased  from  Origene  (-1, 
AGUCGGGAAUUAUAGAUACAGCUTA;  -2, 
CUACUACAAUAAUGACCUGCUCCGA;  -3, 
GGAUAUGCCACAGGACUUCAAGGCT). 


Point  mutation  at  site  of  rsl077858  by  CRISPR 

LNCaP  (AG  genotype  for  rs  1077858)  and  22RV1  (AA  genotype  for  rs  1077858)  cells 
were  chosen  to  reconstitute  rsl077858  GG  genotype  by  CRISPR.  '“Briefly,  5  different 
target  oligos  (Oligo  1-F,  caccgGGC AGGGTCTCGCTCC ATTC ;  Oligo  1-R, 
aaacGAATGGAGCGAGACCCTGCCc;  Oligo  2-F, 
caccgGCAGGGTCTCGCTCCATTCA;  Oligo  2-R, 
aaacTGAATGGAGCGAGACCCTGCc;  Oligo  3-F, 
caccgGAGGTTACGTAGTCCCTGAA;  Oligo  3-R, 
aaacTTC AGGG ACT ACGT AACCTCc ;  Oligo  4-F, 
caccgTTGTATCTCTCACATTAGAC;  Oligo  4-R, 
aaacGTCTA  AT  GT  GAG  AG  AT  AC  A  Ac ;  Oligo  5-F, 
caccgCGAGACCCTGCCTTGTCTGA;  Oligo  5-R, 
aaacTCAGACAAGGCAGGGTCTCGc)  were  synthesized  and  cloned  into 
lentiCRISPRv2  one  vector  system,  respectively.  A  single  strand  DNA  oligo  of  101 
nucleotides  was  synthesized  as  the  donor  template,  which  have  50bp  of  homology  arms 
on  each  side  flanking  the  change  site  for  rsl077858  A—>G.  The  donor  DNA  sequence  is 
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AAAGGCCCCTGTCAGATGGGGGAGACAGTCTCTGAGCTGGGGGAAGCCCCT 


GTCTAATGTGAGAGATACAATCTCTCTGATGGAGGTTACGTAGTCCCTGAAT 
GGAGCGAGACCCTGCCTTGTCTGATGGAAGAGACCTGACCTCTGTGATGGAG 
AAGGCCCCTGTCTGATAAGGGAGACG.  2jal  of  lOpM  donor  oligo  was  used  for 
each  well  of  24- well  plate  to  co-deliver  into  the  cells  with  5  constructs.  After  puromycin 
selection,  cells  derived  from  single  colony  were  assayed  for  point  mutation  (Taqman 
SNP  assay)  and  SLC02B1  mRNA  levels  (RT-PCR). 

We  have  assayed  22  single  colonies  for  LNCaP  cells  and  109  single  colonies  for 
22RV1  cells.  Unfortunately,  we  did  not  detect  any  A  G  mutation  for  LNCaP  cells.  For 
22RV1  cells,  89  single  colonies  remained  AA  genotype,  15  single  colonies  are  AG,  and 
5  of  them  are  GG. 


1.  Ran  FA,  Hsu  PD,  Lin  CY,  et  al:  Double  nicking  by  RNA-guided  CRISPR 
Cas9  for  enhanced  genome  editing  specificity.  Cell  154:1380-9,  2013 

2.  Sanjana  NE,  Shalem  O,  Zhang  F:  Improved  vectors  and  genome-wide 
libraries  for  CRISPR  screening.  Nat  Methods  11:783-4,  2014 
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Supplementary  Table  1.  Patient  Characteristics  and  Outcomes  in  Initial,  Validation 
and  Combined  Cohorts 


Initial  Cohort 

(N=478) 

N  % 

Validation  Cohort 

(N=616) 

N  % 

Combined  Cohorts 

(N=1094) 

N  % 

At  diagnosis 

PSA,  ng/mL 

412 

14 

537 

9.9 

949 

11.6 

Median,  IQR 

(7.0,46.0) 

(5.7,24.8) 

(6.1,31.8) 

Age,  year 

458 

62 

580 

61 

1038 

61 

Median,  IQR 

(56,67) 

(56,67) 

(56,67) 

Clinical  T  stage 

Tx/Unknown 

194 

40.6 

99 

16.1 

293 

26.8 

T1 

122 

25.5 

379 

61.5 

501 

45.8 

T2 

137 

28.7 

110 

17.9 

247 

22.6 

T3-4 

25 

5.2 

28 

4.5 

53 

4.8 

Clinical  N  stage 

Nx/Unknown 

278 

58.2 

349 

56.6 

627 

57.3 

NO 

166 

34.7 

213 

34.6 

379 

34.6 

N1 

34 

7.1 

54 

8.8 

88 

8.0 

Clinical  M  stage 

Mx/Unknown 

217 

45.4 

316 

51.3 

533 

48.7 

MO 

181 

37.9 

208 

33.8 

389 

35.6 

Ml 

80 

16.7 

92 

14.9 

172 

15.7 
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Initial  Cohort 

(N=478) 

N  % 

Validation  Cohort 

(N=616) 

N  % 

Combined  Cohorts 

(N=1094) 

N  % 

Biopsy  Gleason  score 

6  or  less 

77 

16.1 

81 

13.1 

158 

14.4 

7 

153 

32.0 

198 

32.1 

351 

32.1 

8  or  more 

171 

35.8 

274 

44.5 

445 

40.7 

Unknown 

77 

16.1 

63 

10.2 

140 

12.8 

Type  of  local  therapy 

RP  +/-RT 

187 

39.1 

264 

42.9 

451 

41.2 

RT  only/other 

143 

29.9 

201 

32.6 

344 

31.4 

None 

148 

31.0 

151 

24.5 

299 

27.3 

Use  of  ADT  as  part  of  local  therapy 

84 

17.6 

204 

33.1 

288 

26.3 

At  ADT  initiation 

PSA,  ng/mL 

380 

14.9 

586 

10.1 

966 

11.9 

Median,  IQR 

(5.2,68.4) 

(4.1,36.5) 

(4.5,46.8) 

Presence  of  metastases 

274 

57.3 

360 

58.4 

634 

58.0 

Received  Anti- androgen  during  ADT 

333 

69.7 

414 

67.2 

747 

68.3 

Received  intermittent  ADT 

90 

18.8 

221 

35.9 

311 

28.4 

Abbreviations:  ADT,  androgen  deprivation  therapy;  RP,  radical  prostatectomy;  RT, 
radiation  therapy;  M,  metastatic  disease 
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Supplementary  Table  2.  Association  of  SLC02B1  Genotype  with  TTP  on  ADT  in 
Combined  (Initial  plus  Validation)  Cohort  Stratified  by  Status  of  Prior  Hormone  Therapy 


No  prior  hormone 

With  prior  hormone 

N 

Adjusted  HR 

(95%  Cl) 

N 

Adjusted  HR 

(95%  Cl) 

P  value 

(interaction 

rs  12422 149 

AA/AG 

144 

1  (reference) 

53 

1  (reference) 

0.642 

GG 

659 

1.29(1.04,1.60) 

235 

1.44  (0.97,2.13) 

rs  1789693 

AA 

318 

1  (reference) 

109 

1  (reference) 

0.517 

AT 

358 

0.95  (0.79,  1.13) 

128 

1.17  (0.85,  1.59) 

TT 

123 

0.99  (0.78,  1.26) 

46 

1.09  (0.72,  1.67) 

rs  1077858 

AA/AG 

684 

1  (reference) 

255 

1  (reference) 

0.104 

GG 

116 

1.32(1.06,1.65) 

31 

0.87  (0.55,1.37) 
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Supplementary  Fig  1.  Kaplan-Meier  curves  of  time  to  progression  during  ADT  in  the 
validation  cohort,  according  to  SLC02B1  genotypes  (A)  rsl2422149,  (B)  rsl789693,  (C) 
rsl077858. 


(A) 


GG  506 


309 


188 


116 


68 


43 


(B) 


xi 

ra 

xi 

0 


No.  at  risk 


I4  $T 

Months  since  ADT  initiation 


rsl  789693 
_  AA 

-  AT 

-  TT 


AA 

238 

152 

95 

65 

35 

20 

AT 

281 

171 

107 

60 

41 

29 

TT 

90 

57 

39 

24 

14 

7 

AA/AG  530  339  211  130  76  47 

GG  81  41  28  20  14  9 
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Probability 


Supplementary  Fig  2.  Kaplan-Meier  curves  of  time  to  progression  during  ADT  in  all 
cohorts  (initial  plus  validation  cohort)  according  to  rs  12422 149  genotypes  as  stratified  by 
metastatic  status  (A)  MO  versus  (B)  Ml  at  ADT  initiation. 
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Supplementary  Fig  3.  Kaplan-Meier  curves  of  overall  survival  from  ADT  initiation  in 
all  patients  (initial  plus  validation  cohort),  according  to  SLC02B1  genotypes  (A) 
rs  12422 149,  (B)  rsl789693,  (C)  rsl077858. 


GG  894  841  735  585  440  337  261  190  138  105  74 


No.  at  risk 

Years  from  ADT  initiation 

AA 

427 

399 

348 

283 

207 

158 

117 

79 

55 

43 

30 

AT 

486 

461 

405 

327 

237 

195 

149 

117 

84 

62 

48 

TT 

169 

158 

139 

111 

88 

63 

52 

40 

28 

21 

14 

AA/AG  939  880  771  628  466  365  283  210  151  114  81 

GG  147  141  124  97  70  53  36  25  16  12  10 
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Supplementary  Fig  4.  Kaplan-Meier  curves  of  overall  survival  from  ADT  initiation  in 
all  cohorts,  according  to  SLC02B1  genotypes  at  rsl2422149  (A  and  B)  andrsl077858 
(C  and  D),  stratified  by  metastatic  status  at  ADT  initiation. 


(A)  MO 


AA/AG  93  89  81  73  53  45  34  24  17  12  10 

GG  366  359  322  266  226  179  139  100  79  67  44 


(B)  Ml 


AAAG  104  96  85  71  46  37  25  22  12  9  7 

GG  528  482  413  319  214  158  122  90  59  38  30 


(C)  MO  (D)  Ml 


AA/AG  397  387  347  296  241  194  151  107  84  69  45  AA/AG  542  493  424  332  225  171  132  103  67  45  36 

GG  60  59  54  42  37  30  22  17  12  10  9  GG  87  82  70  55  33  23  14  8  4  2  1 
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